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ABSTRACT 


Border protection is a vital national security issue 
for most countries. The U.S. Customs and Border Protection 
(CBP) is responsible for protecting the borders of the U.S. 
from terrorism, human and drug smuggling and illegal 
migration. The U.S. CBP improves manpower, technology and 
infrastructure along the border through various projects. 

In this study, part of the Tucson sector in Arizona is 
modeled in an agent-based model (MANA) to explore the 


effects of using a hand-launched, mini Unmanned Aerial 





Vehicle (miniUAV) along with other assets, such as Border 





Patrol (BP) agents, surveillance towers, the Predator B, 





seismic sensors and communication centers. 

The results from the runs of different scenarios, 
created by a Nearly-Orthogonal Latin Hypercube (NOLH) 
design, are analyzed using comparison tests, linear 
regression, and regression trees. 

As a result, the use of miniUAVs is found to be 
beneficial in capturing the illegal entrants in this 
analysis and thus could potentially provide more secure 
borders. Adequate manpower, in this case BP agents, and a 
reliable communication web to compose a Common Operational 
Picture (COP) emerge as the most important factors 


regarding border protection in this analysis. 
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EXECUTIVE SUMMARY 


Border protection is a vital national security issue 
for most countries. The U.S. Customs and Border Protection 
(CBP) is responsible for protecting the borders of the U.S. 
from terrorism, human and drug smuggling, and illegal 
migration. The U.S. CBP improves manpower, technology and 
infrastructure along the border through various projects. 

The southwest border of the U.S. contains the majority 
of the illegal activities. In this study, part of the 
Tucson sector in Arizona is modeled in MANA, an agent-based 


model, to explor the effects of using a hand-launched, 





mini Unmanned Aerial Vehicle (miniUAV) along with other 
assets, such as Border Patrol (BP) agents, surveillance 
towers, Predator B, seismic sensors and communication 
center. 

The main focus of the Secure Border Initiative (SBI), 


one of the major projects regarding border protection, is 





the integration of border security programs to gain 


effective control of the U.S. borders through substantial 





investments in technology, infrastructure, and enforcement 
personnel. SBiInet Technology Program is one of the two 
programs under SBI for acquiring, developing, integrating, 
and deploying an appropriate mix of surveillance 
technologies, such as cameras, radars, and sensors, and 


command, control, communications, and intelligence (C3T) 








technologies. 
The elements and concept of SBInet are used to develop 


the border security model in MANA. Ground sensors, capable 





of detecting illegal activities, are the first layer of 





defense. Although ground sensors can not differentiate 


XVLL 





between legitimate and illegal activities, they can 





transmit the detection information to the towers. 


Towers have day and night capable video surveillance 





systems that can detect and classify illegal entrants. 
Construction of a surveillance tower takes years in some 
circumstances due to the land acquisition processes. 

The Predator B is the mid-altitude UAV that is 
modeled. With its long endurance, high speed and state of 


the art airborne sensors, the Predator B is a high value 





asset that is currently used by the U.S. CBP to detect, 


classify and track illegal entrants. 





At the tactical level, the Raven, a hand-launched 





miniUAV, is modeled in MANA, which can be used directly by 
BP agents without the need for any coordination. Although 





it contributes BP agents to classify and track the illegal 
entrants, the miniUAV has some limitations such as wind 
speed, short endurance and operational radius. 

The BP agents patrol in their area of responsibility 


to prevent illegal activities. They can receive information 








from all of the sensors and proceed, mounted or dismounted, 





to the area where an illegal activity is classified by 
sensors. Only the BP agents are capable of capturing 


illegal entrants. 





All the information gathered by sensors is transmitted 





to the BP agents via the communication center. This center 
is the essential part of the SBInet concept, the main focus 
of which is to create a COP. 

Once a base model is created in MANA, a Nearly- 
Orthogonal Latin Hypercube (NOLH) design is used to create 


various scenarios to explor th effects of factors on 








different measures of ffectiveness (MOES). MOEs, 


XVLIiL 


determined for analysis, are number of captured immigrants 


or smugglers and number of classified immigrants or 





smugglers. 

After running the model using 257 design points with 
30 replications of each, 7710 data points are acquired for 
each of four MOEs. JMP is used as a statistical tool to 


examine this data through comparison tests, linear 





regression, and regression trees. In order to fit models, 
the 7710 data points are summarized by their means and 


standard deviations, thus collapsing the data set for 





analysis down to 257 data points. 


According to a comparison test, the use of miniUAVs is 





found to be beneficial in capturing the illegal entrants 
and thus provides more secure borders. The mean values of 


the number of captured illegal immigrants are calculated as 








3.04 and 5.49, out of 20, for the scenarios without miniUAV 


and with miniUAV, respectively. Although there is no 





Significant difference for the number of classified illegal 








entrants, between the scenarios with the miniUAV and 











without mini UAV, it is found that using mini UAV decreases 








the time needed to classify all illegal entrants. 

Adequate manpower, in this case BP agents, and a 
reliable communication web to compose Common Operational 
Picture (COP) emerge as the most important factors 
regarding border protection, according to linear regression 


models. Some other important factors are classification 





ranges of towers, the presenc of the Predator B and 





miniUAV, endurance of the mini UAV, the probability of 


classification of the towers and Predator B, and the tower 





height. 
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z. INTRODUCTION 


A. BACKGROUND 


Border protection is a vital security issue for many 
countries. Terrorists; human, drug, contraband, and weapons 
smugglers; and immigrants have been trying to cross the 
borders illegally for decades. These illegal activities 
cause security problems especially for the regions adjacent 


to borders and are threats for the whole of the country. 


1. Overview of the U.S. Borders 


U.S. shares roughly 7,000 miles of land border with 


Canada and Mexico, 1,900 miles of which is with Mexico 








where most of the illegal activities take place. Customs 





and Border Protection (CBP), under the Defense Homeland 
Security (DHS), is responsible for protecting the U.S. 
borders. According to U.S. CBP, they protect the borders of 





the U.S. from terrorism, human and drug smuggling, illegal 
migration, and agricultural pests while simultaneously 
facilitating the flow of legitimate travel and trade. There 
are 327 Ports of Entry along the U.S. and 144 CBP Border 


Patrol Stations within 20 sectors. The proposed budget of 





the CBP for FY2008 to protect the borders was $8.8 billion 
(Ui Sy CBP, 2008)~ 


2. Organization and Key Offices 


CBP has an organizational structure as depicted in 








Figure 1. Every individual office out of 15 under the CBP 


Commissioner is responsible for various types of issues 





regarding U.S. Border protection. Some of these offices 


that are important for this study are described as follows 


in the U.S. CBP’s Performance and Accountability Report: 


Office of Intelligence and Operations Coordination 
(OIOC): Combining targeting and analysis functions from the 
Office of Field Operations (OFO), the Office of Information 
and Technology (OIT), Office of International Trade (OT) 
and the Office of Border Patrol (OBP), OIOC is responsible 


for the entire intelligence cycle. 








Office of Field Operations (OFO): OFO is responsible 
for conducting 327 Point of Entries (POEs) to prevent 
terrorists and terrorist weapons from entering the United 


States. 


U.S. CUSTOMS & BORDER PROTECTION 


POLICY & PLANNING cea eecin 
DEPUTY COMMISSIONER [omar | 
EQUAL EMPLOYMENT 
OPPORTUNITY Executive Secretariat 


INTELLIGENCE & 
BORDER PATROL FIELD OPERATIONS ee: OPERATIONS 
COORDINATION 
INTERNATIONAL 
AFFAIRS AND TRADE CONGRESSIONAL PUBLIC AFFAIRS CHIEF COUNSEL 
RELATIONS AFFAIRS 


TRAINING & HUMAN RESOURCES INFORMATION & 8ECURE BORDER 
DEVELOPMENT MANAGEMENT TECHNOLOGY INITIATIVE 




















Figure 1. U.S. CBP Organization Chart (From: U.S. CBP, 2008) 


Office of Information and Technology (OIT): OLT 
operates a high-performance Information Technology (IT) 
infrastructure, supports tactical communications, 
scientific solutions, and forensic services, and also 
implements and supports CBP’s TY, automation, and 


technology strategies. 


Office of CBP Air and Marine (A&M): Main missions of 


CBP A&M, the world’s largest civilian law enforcement air 





force, are air and marine interdiction, air and marine law 
enforcement, and air and national border domain security. 
CBP A&M also supports DHS missions such as response and 
recovery to natural disasters and terrorism. CBP A&M 


operates from 74 air and marine locations, with more than 





900 federal agents, 270 aircraft, and more than 180 


maritime vessels. 


Office of Border Patrol (OBP): OBP is the main office 





that is responsible for preventing terrorists, weapons of 





terrorism, illegal aliens, drugs, and smugglers from 





entering the United States between the POES. Organization, 





mission, tactics and capabilities of BP are discussed under 





Border Patrol title in detail. 


Office of Secure Border Initiative (SBI): SBI focuses 


on effective integration of border security programs. There 





are two current programs, SBInet Technology Program and SBI 


Tactical Infrastructure Program, under SBI and various 





projects under these programs. The concept of SBI and its 
components are discussed separately in detail (U.S. CBP, 


2008). 


3: Overview of Illegal Activities 


According to the CBP, 178,770 pounds of cocaine, 2,178 


pounds of heroin, 2,471,931 pounds of marijuana and, 2,770 





pounds of methamphetamine were seized along the border of 





the U.S. in 2008 (U.S. CBP, 2008). Drug smugglers can use 





illegal immigrants as a cover or distraction for their 
ground operations. They are also capable of using many 
kinds of vehicles, aircraft, vessels, and even  semi- 


submersibles that have a low radar profile. 


Illegal immigration is an important issue particularly 








along the southwest border with an average of 1 million 


apprehensions per year. According to Nufiez-Neto, this 





number represents approximately 97% of the all illegal 





alien apprehension along the U.S. borders (Nufiez-Neto, 


2008). 
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Figure 2. Southwest Border Apprehensions (From: Nufiez-Neto, 


2008) 


Nufiez-Neto and his co-authors also state that although 
more than 90% of the aliens trying to enter the U.S. 
illegally between POES are Mexican nationals, the number of 
other than Mexicans (OTMs) being apprehended by the BP has 
more than tripled over the three years between 2002 and 
2004. The majority of these apprehensions have come from 


four nations: Honduras, Brazil, El Salvador, and Guatemala 





(Nufmez-Neto et al., 2005). 


Immigrants have to walk long distances in rough 
terrain and weather conditions while crossing the border 
and these difficulties cause many deaths. According to 


Seghetti, launch of the BP’s “Prevention Through 
5 


Deterrence” strategy in 1995 that focused on pushing 


illegal immigration away from population centers’ led 

















illegal immigrants to attempt to cross the border in remote 





desert regions, increasing the number of migrant deaths. 
Rates increased from 1.6 deaths per 10,000 apprehensions in 
FY1999 to 3.7 ain FY2003 which indicates more than 300 
deaths per year (Seghetti et al., 2005). 


Apprehension rates are decreasing with the help of 
numerous projects and operations that improve technology, 
infrastructure and manpower to detect, identify, track and 


apprehend illegal immigrants, smugglers and contraband. 





These projects also enhance collaboration among the _ key 


agencies and develop new laws to deal with legislative 





issues. 


B. RESEARCH QUESTIONS 


This research is guided by the following questions: 

e How might limited assets (such as UAVs, Border 
Patrol agents, ground sensors, and remote video surveillance 
systems) be configured in order to enhance detection and 


classification rates of illegal border crossings? 








e What are the potential differences in illegal 


border crossing detection and classification when using 





miniUAVs with various capabilities? 


Cc. SCOPE OF THE THESIS 


This thesis discusses prevention of illegal border 
crossings without examining its social and economic effects 
on the U.S. To achieve that, a 30-mile-section of U.S.- 
Mexican border, near Nogales, AZ, is evaluated with various 


entities—Border Patrols, sensors, UAVs, smugglers, illegal 
6 


immigrants, etc.—taken into consideration. The impact of 


various assets on detection and classification of illegal 





border crossings is examined by changing the configuration 
of these assets. This analysis gives insight into how to 


configure different assets used in the detection and 





classification process to explore the interactions of 





assets with each other. 


The reason for choosing Tucson can be understood by 


looking at Figure 3, which displays the Southwest 





apprehension rates by sector. San Diego, El Paso and 





McAllen are the sectors with highest apprehension rates in 


the 1990’s. While rates hav decreased in California and 








Texas, with the help of Operation Gate Keeper in San Diego 
and Operation Hold-the-Line in El Paso, there has been an 
increase in Arizona sectors, especially in Tucson. The main 
reason is the success of the operations and projects of BP 


to push illegal immigration away from the populated areas. 





The shifting of immigration routes from more populated 
areas to the sparsely populated areas or from urban to 
rural is called the “Balloon effect” by Ordéfiez in his 


thesis (Orddfiez, 2006). 
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Figure 3. SW Border Apprehensions by Sector and Fiscal Year. 
(From: Nufiez-Neto, 2008) 


The Tucson sector is the most active region for 


illegal activities. As a result, a majority of the 





technological projects under SBiInet Technology Program are 
initiated in the Tucson sector, such as Project 28, a 


finalized proof-of-concept project, deployed along an 





approximately 28-mile section of the Tucson Sector. 


D. LITERATURE REVIEW 


In his thesis, Pulat (2005) develops an optimization 


model by maximizing the probability of escape of an 





infiltrator and then exploring the U.S. BP’s courses of 
action to minimize the maximum achievable probability of 
escape. This model, focused on the U.S - Mexican border 
near Yuma, Arizona, provides an insight to identify 


critical road segments and areas to defend, and also 


evaluates the effects of employing different types of 


assets and strategies on the infiltration patterns. 


Ordéfiez (2006) illustrates the outcome apprehension 
probability of migrants given the implementation of various 
operational strategies by presenting the Arizona-Sonora 


Border (ASB) Model. 


Patrascu (2007) develops a model in Microsoft Excel 


that maximizes the probability of detecting intruders by a 








distributed sensor network subject to a budgetary 


constraint. 


Beeker and Page (2006) discuss a federation approach 
to simulate southern U.S. Border in MANA. Although there is 
no information about the model itself, the paper explains 
how they solve the issue of modeling a long border by 


dividing into smaller parts that overlap each other. 


In his study, Raffetto (2004) analyzes Marine 





Expeditionary Brigade (MEB) commander’s 2015 unmanned 





aerial vehicle (UAV). He uses MANA to explore the impact of 
various capabilities of the UAV on intelligence gathering 
missions. Through modeling and data mining, this study 
provides insights into the importance of various UAV 
characteristics, such as airspeed, endurance, sweep width, 


and sensor capability. 


THIS PAGE INTENTIONALLY LEFT BLANK 
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AB Bet BORDER SECURITY CONCEPT 


There are three legs of the border security concept: 





infrastructure, technology and manpower. Technology and 











infrastructur are developed by two different projects 








under the Secure Border Initiative (SBI) and Border Patrol 





agents are the main manpower to conduct the operations. 
Many types of sensors, radars, day/night cameras both on 
the ground and on airborne platforms form the technological 
barrier system along with a communication web to obtain 


operational pictures. 


A. SECURE BORDER INITIATIVE 


According to CBP website, the main focus of SBI, 
launched in 2005, is the integration of border security 
programs to gain effective control the U.S. borders through 
substantial investments in technology, infrastructure, and 


enforcement personnel. The Tactical Infrastructure Program 





and SBInet Technology Program are the two current programs 
under SBI (U.S. CBP). 
1. SBI Tactical Infrastructure Program 


According to a Government Accountability Office (GAO) 


Report, the SBI program received about $3.6 billion from 





FY2006 to FY2009 and about $2.4 billion of this amount has 





been allocated to complete approximately 670 miles of 





vehicle and pedestrian fencing along the border between th 
United States and Mexico. This program consists of building 


various kinds of pedestrian fences, vehicl fences, all- 





weather patrol roads and permanent lighting (Stana, 2009). 
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In the same report it is stated that most of the 
pedestrian fencing is single layer fencing and in some 
places, a secondary fencing is constructed parallel to 
existing primary fencing for additional operational 
advantage to deter illegal cross-border activities. 
Pedestrian fences are mostly constructed in urban areas 
(Stana, 2009). An example of the infrastructure system is 
illustrated in Figure 4. Sub-projects under SBI Tactical 


Infrastructure Program are listed in Table l. 























Fence projects Pedestrian fence miles | Vehicle fence miles Total miles 

PF 70 81 N/A 81 
PF 225 210 N/A 210 
VF 300 N/A 227 227 
Legacy pedestrian 

fence 67 N/A 67 
Legacy vehicle fence N/A 76 76 
Total 358 303 661 














Source: SBI. 





Note: N/A = not applicable. 





Seventy-eight miles of pedestrian fencing and 57 miles of vehicle fencing were in place before the SBI 
program began. However, since SBI began construction, some miles of fencing have been removed, 


replaced or retrofitted resulting in mileage totals that are different from those we have previously 


reported. 








Table 1. 


Border Fence Projects 
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(From: 


GAO, 2009) 
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Primary Fence 


All-Weather Patrol Road 


12" x 48° Concrete Footng 





Figure 4. Border Infrastructure System in San Diego (From: 
Nufiez-Neto&Kim, 2008) 


2. SBInet Technology Program 


SBInet, launched in 2005 as a part of SBI, replaced 





two former programs, America’s Shield Initiative (ASI) and 
the Integrated Surveillance Intelligence System (ISIS). A 
GAO report describes SBInet as the program under SBI for 
acquiring, developing, integrating, and deploying an 
appropriate mix of surveillance technologies, such as 
cameras, radars, and sensors, and command, control, 
communications, and intelligence (C3I) technologies (GAO, 


2008). 
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Figure 5. Conceptual Depiction of Long Term SBInet Operations 


(From: GAO, 2008) 


In a review by DHS Office of Inspector General it is 


stated that ISIS was formally established in 1998 as a 





technological assistance program to provide continuous 
monitoring of the borders in all weather conditions. The 
main components of ISIS are sensors, Remote Video 
Surveillance Systems (RVSS) and the Intelligent Computer 


Assisted Detection System (ICAD). Sensors send a signal to 





the ICAD whenever detection occurs and then this area is 
investigated by a Law Enforcement Communication Assistant 
(LECA) by using RVSS or by a BP agent. Recognizing the need 
to improve border surveillance and remote assessment and 
monitoring technology, OBP began developing America’s 


Shield Initiative (ASI) in June 2003, as a program to 





integrate surveillance technology, communications, and 
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visualization tools by maintaining and modernizing ISIS 


(Office of Inspector General, 2005). 


According to CBP website, SBInet develops and deploys 





new integrated technology solutions for front-line CBP 








personnel by providing enhanced detection, tracking, 
response, and situational awareness capabilities. An 
essential element of SBInet Technology Program is a 


Command, Control, Communications, and Intelligence (C31) 


Common Operating Picture (COP), which provides real-time 





Situational awareness of the area of responsibility of BP 
agents and allows fewer agents to cover more ground (U.S. 


CBP). 


Current and completed SBInet projects that are 


mentioned in the CBP website include: 


a. Ajo-1 Project 


The Ajo-1l Project, expected to be fielded in 
calendar year 2010, is located south of Ajo, Arizona and is 
comprised of six sensor towers, Six communication towers, 
200 unattended ground sensors, and 20 miles of access road 


improvements and construction. 


b. Command, Control, Communications, and 
Intelligence (C3I) Project 


This on-going C3I project, launched in December 
2007, focuses on obtaining a Common Operating Picture 
(COP). 

c. El Paso Phase I & Phase II Projects 


These projects include a tower system equipped 


with a camera, radar, sensors and communications equipment; 
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new and retrofitted vehicles; and unattended ground 
sensors. Phase I and phase II covers an approximately 74- 
mile area of the U.S./Mexico international border in the 
vicinity of Fort Hancock, Fabens, and Ysleta, Texas and a 
194-mile area from New Mexico (Lordsburg, Deming, and Santa 


Teresa) into El Paso, Texas. 


d. Northern Border Project 


The Northern Border Project includes, but is not 


limited to, 16 Remote Video Surveillance Systems (11 sites 








in the Border Patrol’s Detroit Sector and 5 sites in the 
Border Patrol’s Buffalo Sector), and 3 Mobile Surveillance 


Systems. 


e. Project 28 


As a finalized proof-of-concept project, Project 
28, deployed along an approximate 28-mile section of the 
Tucson Sector, in the vicinity of Sasabe, Arizona, has 
provided some insights on the operational and technical 


challenges of the SBInet program. 


£. TUS-1 & Tucson Sector Projects 


TUS-1, the first of two operational deployments 


of the SBInet technology solutions, includes nine sensor 





towers; eight communication towers, 200 unattended ground 
sensors, and six miles of access road improvements and 
construction over 23 miles of the U.S-Mexico border near 


the Sasabe Port of Entry. Once key elements of the project 








plan are strengthened, the Tucson Sector Project will cover 











all 262 miles of border in this sector. 
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g. Yuma Sector Project 


Having the same components with the other 
projects, Yuma Sector Project covers approximately a 125- 
mile area in Yuma Sector, between the Camp Grip training 
facility and the El Centro Sector line near the Imperial 


Sand Dunes (U.S. CBP). 
3. Key Components of SBInet 


a. Unattended Ground Sensors (UGS): 


Seismic and magnetic UGSs construct the primary 


layer for detection of illegal activities. There are mainly 





two types of ground sensors; older technology sensors, with 


high false alarm rates, that are deployed as a part of ISIS 





and new technology sensors as a part of SBInet. 


For older sensors a report by the Office of 





Inspector General under DHS states that, these sensors 
are the most used as well as the easiest and least 
expensive to install and maintain. The sensor sensitivity 


level can be adjusted to help filter false alerts. There 





are more than 11,000 sensors along the northern and 





southwest borders (Office of Inspector General, 2005). 


It is stated in a hearing before the U.S. Senate 





that averag detection range for ground sensors is 10 


meters (Senate Hearings, 2005). On the other hand, Dumpert 








and Dirksen claim 30 feet for the detection range. (Dumpert 
& Dirksen, 2006) Although effective in detecting activity 


or movement, these sensors cannot differentiate between 








illegal activity and legitimate events. Therefore, very 
detection needs to be investigated by either a BP agent or 


some other sensor that has classification capability. But, 








due to high false alarm rates caused by other than illegal 
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alien activity, such as local traffic, outbound traffic, a 
train, or animals, this effort, to clarify detections by 
UGSs, is mostly ineffective. According to the analysis of 
sensor alerts along both U.S. borders in the same report by 


the Office of Inspector General, 90 percent or more of 





these alerts were false alarms. Data conducted by the 








analyst to determine this result is shown in Table 2. 


Number of Staging, Identified | Unidentified, 
Ticket ICAD Turn or False Unknown or Not 
Source Tickets Apprehensions | Got Away | Alarm Available 
Sensor 
Alerts 29,710 252 % 








RVS Camera 
Observations 





Non-ISIS 
Sources 780 382 49% 
Source: OIG analysis of OBP ICAD report data. 

Note: Rows may not equal 100 percent due to rounding. 


























Table 2. Southwest Border ICAD Ticket Results (From: 
Office of Inspector General, 2005) 


Operation of a UGS is affected by moisture which 
can cause corrosion, and intentional or accidental physical 
damage by vehicles, machinery, or vandals. Also insects 
penetrating sensors can affect functionality (Office of 
Inspector General, 2005). 

According to Knobler & Winston, due to these 
issues of the legacy sensors, a derivative of OmniSense® 
sensor system, developed by McQ, with terrestrial 
communications and no imaging capabilities has been 
delivered as part of the SBInet program. Five modalities of 
transducers are used in CORE sensors are seismic, acoustic, 


Passive infrared (PIR), magnetic, and tripwire. McQ 
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developed a fusion architecture and protocol for passing 
the appropriate data among devices to achieve the following 
goals: 

° Reduction of false alarms compared with 
detections from individual sensors 


* Improved classification compared with detections 





from individual sensors 
* The ability to track target speed and direction 
(Knobler & Winston, 2008). 


b. Remote Video Surveillance Systems (RVSS) 


According to a GAO report, while radars mounted 
on fixed and mobile towers detects movement, cameras on 


fixed and mobile towers are used to identify, classify, and 








track items of interest detected by the ground sensors and 





the radars (GAO, 2008). 


According to the U.S. CBP, these systems are 





deployed on diverse platforms such as Static Remote Video 





Systems that are located on towers, buildings, and in some 
areas, poles, Mobile Remote Video Systems (MRVS) that are 
located almost exclusively on scope trucks with the 
exception of the tripod types and the man portable systems, 


and Trailer Remote Video Systems (TRVS) that are generally 





mounted on “Sky Watch” trailer systems (U.S. CBP, 2008). 


Some characteristics of Remote Video Surveillance 





systems (RVSS) are specified by the Office of Inspector 





General (2005, p. 4) as follows: 


RVS systems provide the primary remote 

identification capability. RVS components include 

cameras, mounting poles, radio, and equipment, 

such as cabling and equipment enclosures. The RVS 

system includes both color (day) and thermal- 
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infrared (night) cameras, which are mounted on 
sixty or eighty-foot poles or other structures. 
RVS camera signals are transmitted to the OBP 
sector or station communications center via a 
wireless system such as microwave signal, or, in 
one sector, via fiber optic cable. Personnel at 
designated communications centers can control 
most RVS cameras remotely using toggling 
keyboards. There are 255 operational RVS camera 
sites along the northern and southwest borders. 











David Aguilar, Chief, Office of Border Patrol 





(OBP), stated in a hearing before the U.S. Senate that out 
of these 255 camera sites, Tucson Sector has 39, Yuma has 


18, Swanton has six, El Centro Sector has 41, El Paso has 





27, Laredo has 20 and, McAllen has 29 and each one of these 


poles has the capability of looking in either direction 





about six miles (Senate Hearings, 2005). 





One of the problems with stable towers, he 





mentions in his speech, is that approximately 92 percent of 





the border is environmentally sensitive in Arizona. 
Therefore it takes a multi-year process to even plant a 
pole in the ground for an RVS camera, to build the tactical 
infrastructure or to build the roadways. Also some towers 
must be constructed in places that are feasible in terms of 


acquisition process instead of the places that provide the 





best coverage (Senate Hearings, 2005). 





Another issue stated by the Office of Inspector 
General is that RVS cameras do not have the ability to 
detect movement by themselves. It should be monitored by a 
person, called a Law Enforcement Communication Assistant 
(LECA), in a BP station or in a command center. In case of 


a triggered alarm by a UGS, the LECA is required to select 





the appropriate RVS camera, manually maneuver the camera in 
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the direction of the sensor and attempt to identify the 
cause of the sensor alert. However, the number of LECAs is 
inadequate. For instance, at one location they visited, 
only one LECA was on duty performing radio-dispatch duties, 
processing sensor alert information via ICAD, and 


monitoring 32 cameras (Office of Inspector General, 2005). 


McDaniel states the specific requirements of 


EO/IR sensors as follows: 


e Detection of Man-Sized Targets at 3500 


meters (minimum) 


e Recognition of Man-Sized Targets at 2200 


meters (minimum) 





e Detection of Vehicle-Sized Targets at 7500 





meters (minimum) (McDaniel, Hughes, & 
Seibel, 2006). 





Figure 6. Tower Deployed in Tucson Sector with Camera and 


Radar (From: GAO, 2008) 


wal 


(oe Unmanned Aerial Vehicle 


Another asset that supports Border operations is 
Unmanned Aerial Vehicles (UAVS) that are capable of 


detecting, classifying and tracking illegal activities. 





According to CBP, CBP A&M currently operates five Predator 
B UAVs; four are from the UAS Operations Center in Sierra 


Vista, Arizona, and one is from the UAS Operations Center 





in Grand Forks, North Dakota (U.S. CBP). The history of the 


UAV demonstrations sponsored by DHS is displayed in Table3. 
























































Unmanned 
Aircraft Sponsor Sorties | Hours 
Demonstration | Location Used (Support) Dates Flown Flown 
Operation Gila ICE (Air Oct-Nov 
Safeguard Bend, AZ PredatorB | Force) 2003 15 106 
King 
Salmon, USCG Nov 
Alaska Demo 1 | AK Predator (Navy) 2003 5 iS £0) 
King 
Salmon, USCG Aug 
Alaska Demo 2 | AK Altair (NASA) 2004 3 36 
Wallops 
Island, USCG Nov-Dec 
VA Aerosonde (NASA) 2004 
Sierra Hermes CBP Jun-Sep 
ABCI Vista, AZ 450 (Navy) 2004 65 590.1 
Nov 
2004 -— 
ABCI Follow- Sierra CBP Jan 
on Vista, AZ Hunter (Army) 2005 4l 329.1 
Borinquen Feb 
Coastal Areas |, PR Aerosonde | USCG 2005 
Table 3. DHS-Sponsored Unmanned Aircraft Demonstrations 
(From Unmanned Systems Roadmap 2007-2032, 2007) 
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According to Bolkcom and Nufniez-Neto, some 
benefits of using UAVs in border operations are that they 


could fill a gap in current border surveillance by 





improving coverage along remote sections of the U.S. 
border, UAVS can provide precise and real-time imagery to a 
ground control operator. UAVs also would have a greater 


chance of tracking an illegal entrant in dense woods or 





mountainous terrain, with thermal detection sensors than 
the stationary video equipment which is often used on the 


borders. 


Despite potential benefits of using UAVs for 


homeland security, there are some common issues such as the 





high accident rate of UAVs, which is currently 100 times 
higher than that of manned aircraft, and UAVs equipped with 
only an EO and IR camera can suffer from cloudy conditions 
and high humidity climates that can distort the imagery 
produced by EO and IR equipment (Bolkcom & Nufiez-Neto, 
2008). 


(1) Predator  B: According to CBP, the 





Predator B (Figure 7) is a tactical multi-mission UAV and 





the primary mission equipment consists of a Raytheon 
AN/AAS-52(V) MTS EO/IR sensor turret/laser designator and a 
General Atomics AN/APY-8 Lynx SAR (U.S. CBP, 2009). 


23 








Figure 7. Predator B (From: U.S. CBP) 

















Max level speed 240 kt (444 km/h; 276 mph) 

Cruising speed 210 kt (389 km/h; 242 mph) 

Max operating 15,545 m (51,000 ft) 

altitude 

Mission radius 1,655 n miles (3,065 km; 1,904 
miles) 

Max range 4,950 n miles (9,167 km; 5,696 
miles) 

Endurance 35 h 

Table 4. Predator B Performance(From: Jane’s) 


The AN/AAS-52 Multi-Spectral Targeting 
System (MTS), one of the mission components of Predator B, 
combines electro-optical (EO), infrared (IR), laser 
designation, and laser illumination capabilities in a 
Single sensor package. The multi-use system offers long- 
range surveillance, target acquisition and tracking 


(Raytheon Company). 





Parameters Features Parameters Features 

60 degrees up, -120 
Wide: 34-45 Elevation: degrees down 

Gimbal slew 
edium-wide: 17 x 22 | rate: 3 radians/sec elevation 
aximum air 
edium: 5.7 x 7.6 speed: >350kts IAS 
edium-narrow: 2.8 x 
Spo Image Fusion: Included 

arrow: 1.2 x l. Automatic video | Multimode (centroid, 
IR and TV) tracker: area, and feature) 
Ultra-narrow: 0. Compliant with MIL-E 
0.8 (IR) Environmental: 5400, MIL-STD-810 
Ultra-narrow: , 1553 data bus and/or 
0.27 CLV) Interface: discrete controls 

Jeol Cie al 0 BG i> ee Le RS-170 (525-line), 

(IR), digital, 

0.11 x 0.14 (TV) Video outputs: other formats available 
Aly OES CX 0.2 
(IR), 
0.06 x 0.07 (TV) Cooling: Self contained 

ultiple sensors such as 
EO-TV, 

image intensified TV, 
Azimuth: 360 illuminator, eyesafe 
degrees, continuous rangefinder, 

spot tracker, image 
fusion, 

and other avionics 


Table 5. AN/AAS-52 MTS Specifications (From: Raytheon 
Company Web Site) 




















Fields of 
view, 
d Ss: 





























Electronic 
zoom, IR & 
TV: 














Gimbal 
angular 
coverage: 














Another component of Predator B is the high- 
resolution Lynx Synthetic Aperture Radar (SAR). In their 





study, Tsunoda, et al., state that the Lynx SAR weighs less 





than 120 lbs., has a slant range of 30 km, allows two 


resolution selections and is capable of 0.1 m resolution in 





spotlight mode and 0.3 m resolution in stripmap mode. In 





ground moving target indicator mode, the minimum detectable 





velocity is six knots with a minimum target cross-section 








of 10 dBsm. The Lynx user interface features a view manager 
that allows it to pan and zoom like a video camera. 


(Tsunoda et al., 1999) Detailed specifications of stripmap 





and spotlight mode are demonstrated in Table 6 and Table 7. 
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Resolution 0.3 to 3.0 m Both slant range and azimuth 
aaage 7 to 30 km Slant range (3-60 km at reduced 
performance) 
eee 2600 pigeis Only with 16-node system (to 3500 
pixels at coarser resolutions) 
wie Sige 630° At O..3 m resolution, 45 deg. 
depression 
ef=. (45 to 185) equity 18 differ me b sau n sc n 
Squint angle Ae center-line and aircraft velocity 
3 vector 
Table 6. Stripmap SAR Mode Specifications (From: Hensley 
et al., 1999) 
Resolution 0.1 to 3.0 m Select one of five 
Boviee 4 to 25 km Slant range (3-60 km at 
reduced performance) 
Patch Size 2 x (640 x 480) pixels 
View size 640 x 480 pixels Over NTSC video link 
+/- (50 to 130) deg 
Squint angle 5 
a=: (ie te 445) aes 0.15 m resolution and 
coarser 
Table 7. Spotlight SAR Mode Specifications (From: Hensley 


et al., 1999) 


(2) Raven: Army Unmanned Aircraft System 
Operations Field Manual describes the Raven as a man- 
portable, hand-launched small unmanned aerial vehicle 


(SUAV) system designed for Reconnaissance and Surveillance 


(R&S) and remote monitoring. 





Figure 8. Raven B (left), An Operator and Observer Team 
Assembling the SUAS (right). (From: Lifschitz et al., 2007) 
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Raven consists of three main components: UAV 


and sensors, Remote Video Terminal (RVT) and Ground Control 





Unit (GCU). For sensors it has three different 


interchangeable noses, one with front and side look CCD 





color videos and two for infrared cameras, one for front 
other for side look. Specifications of these cameras can be 
found in Table 9. It can easily be assembled in less than 3 
minutes and reusable for more than 100 flights 
(Headquarters, Department of the Army, 2006). 

There are many sources that give slightly 
different numbers for specifications of the Raven but 


according to the Directory of U.S. Military Rockets and 





Missiles it weighs about 1.9 kg (4.2 1b), has a flight 





endurance of 80 minutes and an effective operational radius 
of about 10 km (6.2 miles). Flying speed is 45-95 km/h (28- 
60 mph) at typical operating altitude between 30 and 300 





meters (100-1000 ft). The Raven can be ither remotely 








controlled from the ground station or fly completely 
autonomous missions using GPS waypoint navigation (Parsch, 


2004). 
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UAV Characteristics 


Feature Design Specification Feature Design Specification 
Li-Ion 

Power rechargeable Airspeed 
battery 





23 kt loiter, 34 kt 
cruise, 60 kt dash 





150-1,000 ft (45.72- 


Wing Span Ae). BES Cl. 377 * m1) Altitude 304.8 m) AGL 





4 1b (1.81 kg) 60 to 90 minutes (Li-Ion 


(i215 a4 ka || Meurenee rechargeable) 
with carrying 
case) Payload(s) EO/IR sensors 














Hand-launched/auto land 
17 lb (7.71 kg) Launch/Recovery | recovery on soft, 
unimproved surface 














Two MOS nonspecific 




















8-12 km Crew Sa RAes 
Sensor Characteristics 
Feature Design EO IR 
Pixels 768H X 494V 160H X 120V 
Payload Nose 
mevgne 6.2 oz 6.5 oz 
Table 8. Raven UAV and Sensor Characteristics (From: U.S. 





Army, 2006) 


An OEF/OIF study of close combat missions, 
using small unmanned aircraft systems presents some 
limitations of a small UAV, such as Raven, as follows. 


Because of its small size and required low 





overall weight, Raven has limitations in the imaging 
payloads it can carry. The best arrangement for cameras to 


view images at given slant ranges result in a narrow field 





of view. The test results show that to recognize man-sized 





objects, it should fly at 300 feet AGL or lower. 

Another issue is the lack of stabilizing 
image capability. Even the slightest wind shift causes the 
vehicle to move in all three axes which produce images that 
when viewed through the video screen appear as a bobbing 


effect with the viewed image moving constantly. 
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The information provided by Raven is of 
limited use due to its short endurance, small area of 
coverage and inability to communicate with other 


information systems. 





During the operational test, 21 of 61 hand 


launch attempts ended in failure and nineteen of the failed 





launch attempts occurred during nighttime and in zero to 


calm wind conditions (Lifschitz et al., 2007). 





Despite all the limitations stated above, a 
small UAV like Raven is still promising in border 


operations with its inexpensive cost, lower coordination 





and personnel requirements in comparison with Predator B. 








It can be used by every individual mission package of BP 





agents, at the tactical level, as a forward eye with its 





airborne sensors, for detection, classification and 





tracking of illegal entrants. 


B. BORDER PATROL 


According to the CBP, the United States Border Patrol 


(BP) is the mobile, uniformed law enforcement arm of the 





U.S. CBP within the DHS. The BP, with more than 17,000 


agents in FY2008, is charged with securing the U.S. 











international land border between POE, detecting and 
preventing the entry of terrorists, weapons of mass 
destruction, and illegal aliens into the country, and 








interdicting drug smugglers and other criminals along the 





border. The BP achieves this goal by maintaining 
surveillance, following up leads, responding to electronic 
sensor alarms and aircraft sightings, and interpreting and 


following tracks. 
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In some regions, due to diverse terrain features such 
as uninhabited deserts, canyons, or mountains, the BP 


utilizes a variety of equipment and methods to accomplish 








its mission. Electronic sensors, video monitors and night 
vision scopes are all used to detect illegal entries and 
agents patrol the border in vehicles, boats, manned and 
unmanned aircraft, and afoot. In some areas, the Border 


Patrol even employs horses, all-terrain motorcycles, 





bicycles, and snowmobiles (U.S. CBP, 2009). 


BP conducts its mission within 20 sectors with 140 
Border Patrol stations nationwide and 34 permanent 


checkpoints. 


According to a GAO report, on the southwest border, 
the Border Patrol uses an integrated, multilayered border 
enforcement strategy. Along the border, between official 
ports of entry, are the first two layers, consisting of a 


first called line watch and a second, called line patrol. 





Together, these are where the majority of the U.S. Border 





Patrol agents are deployed, to maintain a high profile to 





deter, turn back, or arrest anyone attempting to illegally 





enter the country. The line patrol layer consists of 


smaller contingents of agents deployed behind the line 





watch units to provide direct support of the line watch 
units. The Border Patrol allocates personnel based on a 


combination of intelligence information about potential 





threats from terrorists and contraband smugglers, as well 





as on the estimated volume of illegal entries. In addition, 





a third layer of enforcement is composed of interior 
traffic checkpoints at which Border Patrol agents monitor 


and stop vehicles at  checkpoints-—both permanent and 
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tactical (temporary)—on major U.S. highways and secondary 
roads that are generally 25 to 75 miles inland from the 


border (GAO, 2005). 


Due to the scope of this study BP structure and 


capability in Tucson, AZ region is explained in detail. 
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Figure 9. Tucson Sector Map (From: CBP) 
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The CBP website states that the Tucson Sector is the 
busiest sector on the Southwest Border with more than 2,900 


agents and covering 262 miles of linear border from the 





Yuma county line to the Arizona/New Mexico state line. The 





eight stations, that the sector contains, are located in 
Ajo, Casa Grande, Tucson, Nogales, Wilcox, Sonoita, Naco, 


and Douglas (Figure 9) (U.S. CBP, 2008). 


Apprehension rates currently are high in Tucson while 
rates have decreased in the other problematic sectors such 
as San Diego, CA and El Paso, TX. David Aguilar, Chief, 
Office of Border Patrol (OBP) has noted the shift of the 


illegal routes from urban to rural. This shift gives 
eal 


criminal organizations a wider array of action and also 





forces BP agents to deal with the rural environment, rural 
dynamics and, much broader area of operations with same 
personnel, technology, and tactical infrastructure (Senate 


Hearings, 2005). 


Aguilar also mentions two typical characteristics of 


Arizona, unlike California and Texas, that hinder control 





of the border. One of them is the absence of checkpoints 





which are a very ffectiv way of controlling major 
roadways. In Arizona, they are forced to use temporary or 


mobile checkpoints that are not as effective as permanent 





checkpoints due to the fact they do not have permanent 
infrastructure such as staging areas and detention centers. 
The second reason for an inefficient border control in 
Arizona is the environmental constraint due to the Federal 


ownership of the land. Any action to construct an 





infrastructure becomes a major Federal action subject and 
takes more than years to implement. The environmental 


sensitivity of the land is another issue. While the BP 





agents has to operate carefully to prevent damage to these 


sensitive areas, by not using cars for instance, the 





smugglers or immigrants can easily use these areas without 


any concern. This creates a disadvantage to the BP agents 





in tracking illegal entrants (Senate Hearings, 2005). 
Cc. ILLEGAL ACTIVITIES 
Illegal immigrants usually cross the border with the 





help of human smugglers, in other words, “Coyotes”. Orddéfiez 





states that illegal immigrants meet coyotes in staging 





areas before crossing the border. Coyotes decide which 


gateway to use based upon information provided by 
22 


personnel, known as “Scouts”, at observation posts in the 
border region to provide warnings to the Coyotes (Orddéfiez, 


2006). 


According to Ordéfiez, smugglers and illegal immigrants 





use three main corridors in the BP Tucson Sector to cross 


the U.S. borders: the West Corridor, Nogales Corridor and 





the Naco-Douglas Corridor. These main corridors contain the 





13 most frequently used routes as illustrated in Figure 10: 


* West Corridor: Organ Pipe National Park, Tohono 
O’odham Nation - San Miguel Port of Entry and the Tohono 


O’odham Nation - Pozos Verde. 


* Nogales Corridor: Coronado National Forest, Mariposa 


Canyon, the Nogales POEs and Patagonia. 


* Naco-Douglas Corridor: Fort Huachuca Mountains or 
Copper Canyon, the Naco POE, San Pedro Valley, Whitewater 
Wash area, the Douglas POE and the San Bernardino Valley 
(Ordéfiez, 2006). 
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Figure 10. Identified Most Frequently Used Illegal Human and 
Drug Smuggling Routes (From: Orddfiez, 2006) 


Ordéfiez also states that drug and human smugglers 
usually respect their area of operation and use different 
routes. Nevertheless, in some cases they intersect with 
each other and use this as an advantage. For instance, drug 


smugglers use illegal immigrants to distract the BP agents 





(Ordofiez, 2006). Table 9 shows the classification of 13 





smuggling routes as drug and human. 
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ROUTES 
ROUTES (from E to W along border) MIGRANT 
West Corridor (L3) 
Organ Pipe National Park 
TO Nation —- San Miguel 
TO Nation - Pozos Verde 
Nogales Corridor 
Coronado National Forest (Arivaca) (L3) 
Mariposa Canyon (L2) 
Nogales POE (L1) 
Patagonia (L2) 









































Naco-Douglas Corridor 
Ft. Huachuca Mountains (Copper Canyon) 
(L3) 
Naco POE (L2) 
San Pedro Valley (L3) 
White Water (L2) 
Douglas POE (L1) 


























San Bernardino Valley (L2) 





Table 9. Identification and Classification of Illegal 
Routes (From: Orddfiez, 2006) 








35 


THIS PAGE INTENTIONALLY LEFT BLANK 


36 


III. MODEL 


A. WHY MANA 


MANA (Map Aware Non-Uniform Automata) version 4.04.1, 





an agent-based model, is used to simulate the border in 
this study. With the help of emerging behaviors of agents, 


agent-based models are beneficial to simulate complex 





systems, which in particular contain humans as elements. 


Once general characteristics and propensities of the agents 





ar defined, ach agent decides how to act through the 


Simulation to reach their goals. MANA allows us to create a 





border security model with BP agents, immigrants and 


smugglers. 


Another reason for using MANA is the built-in feature 


for communication among agents. Due to the main focus of 





SBInet Technology Program as the integration of all these 


sensors to obtain Common Operation Picture (COP), 





communication is a critical element of Border Security 
Model, and communication links can be easily constructed in 


MANA by this built-in feature. 


Once elevation and terrain maps are loaded, MANA 


calculates the line of sight (LOS) for agents automatically 








which is an important issue to model surveillance systems. 


Finally, MANA as a distillation model is easy to 
understand for a first-time user and creating a scenario 


from scratch does not take too much time. 


However, besides all these advantages, MANA has’ some 
limitations on simulating some BP operations that are 


different from military applications. Also it is not 
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supported with sufficient analysis tools which cause an 


extra effort to interpret results. 


B. MANA IN GENERAL 


This part is a brief summary of the MANA user manual. 





Refer to the manual for detailed description and features 


of the model. 


The main elements of the MANA model are squads and 
agents that construct these squads. A squad consists of 
either an individual agent or a group of agents with the 


same characteristics. Each squad can be either friend, 





enemy or neutral. 


Each squad has some characteristics associated with 


its agents such as personality and physical capabilities 





and also characteristics of the assets that the agents use 


such as sensors, weapons and communication tools. 


The agents move in a battlefield consisting of a two- 


dimensional grid. The movement algorithm is a critical 





feature that constructs the emerging behaviors by using 








some penalty calculations. The movement constraints and the 
propensities of the agents to reach some goals such as 
going through the next waypoint, enemy, friend, unknown, 
etc., have the main effect on the direction of this 


movement. “Going” feature of the terrain map and speed of 





the agent determines the velocity. The agents that will be 








considered in that calculation ar determined by 





Situational awareness maps of the sensors of the agent 
organically or via communication link, inorganically. 
Characteristics and propensities of an agent can be altered 


through the simulation by using “trigger states.” For 
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instance, a state for “fuel out” can be defined and speed 
of the agent can be set to zero for this state. In another 
example, a state can be defined for “enemy contact” and the 
propensity of the agent for going through the enemy can be 


increased while that for following waypoints is decreased. 


The sensor of an agent has the capability of detection 


and classification within some ranges. There are simple and 





advance modes for the sensors. In simple mode, a cookie- 


cutter sensor is used for detection and classification. In 








advance mode, the user can define a cone-shaped field of 
view (FOV) with an associated slew rate and also define 
probabilities in a decreasing order by range. The first 


requirement for the detection is to have LOS between two 





agents. An elevation map and the concealment values of the 
terrain and the agent itself are the critical factors in 


LOS calculations. 


The agents can be given weapons to kill each other 
within some range and with a probability of a hit. This 
feature is only used to represent capturing of illegal 


entrants by BP agents. 





There two types of communication links, organic and 





inorganic. Organic or intra-squad communication links are 


between the agents of the same squad and have relatively 





basic settings. Inorganic or inter-squad communication 


links are between the squads and have more options to set 








such as range, latency, capacity, reliability, accuracy, 


etc. (McIntosh et al., 2007). 
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Cc. BATTLEFIELD 








In MANA, the Battlefield is the surface where agents 
perform their missions and they cannot jump out of its 


boundaries. The Battlefield is divided by a _ predefined 





grid, each cell of which can be occupied by a single agent. 
For this study, a 1000x1000 grid is used as the battlefield 
to represent the 50x50 km border in the the Patagonia 
region of the Tucson sector, and every cell represents 


50x50 m in the real world. 


Three types of maps are used in MANA to form the 
battlefield: Background, Terrain and Elevation maps. 
Background is the user interface for the battlefield that 
has no effect on calculations. Terrain and Elevation map 
affect both speed and Line of Sight (LOS) calculations of 
the agents and corresponding sensors. Each type of terrain, 
represented with different RGB color codes, has specific 
characteristics in terms of going, cover and concealment. 
“Going” effects agent’s speed, “cover” effects shot 


probability and “concealment” effect sensor capability for 








detection and classification. Terrain types and 


characteristics can be found in Figure 11. 


Going Cover Conceal Red Green Blue 
1.00 0.00 0.00 0 0 0 
0.00 “1.00 1.00 192 192 192 
0.90 0.10 0.95 64 64 64 


i 
[J 
1.00 0.00 0.00 255 255 0 7 
LighBush | 0.75 0.10 0.30 10 255 10 = 
[DenseBush |7.00 1.00 1.00 40 180 40 | 
Mourtain 0.20 0.00 0.75 59 23 18 fal 


0.75 0.00 0.50 200 200 160 





Figure 11. Terrain Properties 
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The elevation map consists of black, white and shades 
of grey. The highest point is represented as white, while 
black is used for lowest point and shades of grey are used 


in between. For the battlefield of this study “Real World 





Elevation Range Max” is set to 1153 m, the difference 





between the highest and lowest point, where the highest 
point is 2588 m and the lowest point is 1435 m. This 
elevation range is converted to a range of 0 to 225 by MANA 


behind the scenes for LOS calculations. The terrain and 





elevation maps used for the model are presented in Figure 


12. 





(31.25,-110.35) 





Figure 12. Terrain and Elevation Map (Created From: National 
Geospatial Intelligence Agency DTED) 








The Elevation Map is constructed by viewing the 


Digital Terrain Elevation Data (DTED) level 2, from the 











National Geospatial Intelligence Agency, using Falcon View 








v3.3. After noting the difference between the highest and 





lowest points of elevation, the image is converted into a 


grey scale bitmap that MANA would understand. 
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While constructing the Terrain Map, the Naval Research 
Lab’s Geospatial Information Data Base (GIDB) tool v2.2 is 


used to view the roads and other major terrain features of 





the area. Then the layers of available data are converted 
into a bitmap that MANA would understand using the open- 
source image editor, GIMP. As a final step, due to having 


only the information on the major hard-surface roads, the 





roads are extrapolated by adding artificial roads in the 


terrain file. 





The time step, which is three seconds, is not defined 











explicitly. The model, however, it determines all speed, 
rate and run time values implicitly. In accordance with the 
three second-time step, the run time for the model is 








10,800 time steps that corresponds to nine hours in the 





real world. This number is regarded as an adequate time for 





the illegal entrants to reach a secure region in the U.S. 


Excluding some of the assets, such as mobile 
surveillance systems and U.S. A&M support, the SBInet 


concept is reflected in the model as shown in Figure 13. 





The battlefield is divided into responsibility areas for 
each BP car. Ground sensors are placed as the first layer 
of defense in locations where they can work in accordance 


with RVS Towers. Communication Centers are located inside 





the BP AOR Lo recelv detection and classification 





information from sensors and transmit that information to 
the BP agents. The Predator B is modeled as an airborne 
sensor that is responsible for a larger area than the 
battlefield. The MiniUAV, Raven, is also modeled as an 
airborne sensor that is used by an imagery BP agent in each 


BP car. Scouts are located along the border to provide 


42 


information to the smugglers. The illegal immigrants and 





the smugglers are spread along the Mexican side of the 
border heading north. Finally, neutrals spread around the 


Battlefield, crossing the border back and forth, to cause 





false alarms. Figure 13 depicts the locations of the 


agents. 





Figure 13. Reflection of SBInet Concept into MANA Model 
D. AGENTS 
1. Command and Communication Center 


Command and Communication Centers are responsible for 
receiving and processing all information from the remote 
sensors and then transmitting that information to the BP 
agents or more appropriate sensors. There are two CCC’s, 


each to cover half of the battlefield with 27.5 km 





communication range, to ensure that BP agents stay in their 
Area of Responsibility (AOR). Three seconds is designated 
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for the communication delay for the base model; however, a 


broader range of values is used in experiments involving 


important factors such as latency, reliability and 
accuracy. 
2. Remote Video Surveillance (RVS) Towers 


RVS Towers are capable of detecting, classifying and 


tracking illegal activities within their sensor range. A 








RVS Tower can explore the type of activity by receiving 
detection information from ground sensors. Sensor of a RVS 
Tower, located 80 m above ground level (AGL) has a cone 
shaped Field of View (FOV) and can observe 360° around 
itself by some slew rate. It also is capable of zooming in 
and out to classify any detection with the help of trigger 


states. 


There are two RVS Tower packages, four towers in each, 


one of which is responsible for the east half of the border 





and the other for the west. They can communicate with a 








corresponding Communication Center in their region with a 





three-second communication delay in the base model. 


3. Predator B 


There is one Predator B that is capable of detecting, 


classifying and tracking illegal entrants. It can feed 





information to BP via Communication Center. Predator B has 
a larger AOR than 50 km, which the model covers. Therefore, 
it ais given a route that consists of zigzag shaped 


waypoints to represent flying in adjacent regions out of 





the battlefield. Predator B follows the waypoints 





initially. Once it detects an activity either by its own 
sensor or by other sensors, it then proceeds to the 
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detection area to classify and track the illegal entrants. 





It can be refueled by a MiniUAV within 500 m to jump into a 





state that forces it to leave the current location and go 


for another mission. This prevents the use of two critical 





assets, the MiniUAV and Predator B, for the same target. 


4. Unattended Ground Sensor (UGS) Agent 


There are five UGS packages, each of which consists of 





50 individual ground sensors. Having a 10-meter-detection 
range for every individual ground sensor, each UGS package 
covers one kilometer of the border. The UGS agent is given 


a very small cone shaped FOV to achieve this one kilometer 








detection line, and it only has the ability to detect. They 
need to co-operate with the RVS Tower to classify the 
detection. The UGS agents are invisible to the other 


sensors. There are 19 possible locations for UGS and five 








out of these 19 locations are selected randomly for each 
run. 
5. Border Patrol Car 


BP agents, “embussed” by BP car agent, patrol their 


AOR by driving their cars at the beginning. The term 





“embuss” is used in MANA to create parent-child relation. 





Child squads move together with their parent squad until 





they are released. 


There are two BP car agents in the model, one for each 





region. Initially, they patrol their AOR by following their 








way points in the default state. Onc they receiv a 





detection of an illegal activity by their inorganic 


communication web, they fall into one of the contact states 





and proceed to the activity area. Due to inclination of 
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driving off-road while going for the enemy, the BP Car 


agent is assigned a lower speed than the actual. After 








coming closer to the area of interest or detecting the 





entrant by its own sensor, the BP Car releases the BP agent 








and falls into an inactive, terminating state. 


6. MiniUAV (Raven) 


The MiniUAV agent, as an aircraft, has an airborne 
sensor that is capable of observing 360° around. It can 


detect, classify and track any activity within its sensor 





range. Initially the MiniUAV is embussed by the BP agent. 


Once the BP agent is dismounted near the area of interest 











and loses contact with the illegal entrant, it releases th 


MiniUAV to re-establish the contact. 





Due to the high speed advantage and airborne sensor, 
150 m AGL, the MiniUAV can detect, classify and track 


illegal entrants and provide direct, instantaneous image of 





the area of interest within a 10 km communication range. 
The MiniUAV goes toward the enemy as long as it has a 


contact, whether organic or inorganic. It also has a patrol 





route to search through in case of no contact. The MiniUAV 
flies one hour before it runs out of fuel and becomes 


permanently inactive. There is one MiniUAV for each BP car. 


Ws Border Patrol (BP) Agent 


Initially embussed by the BP Car, the BP agent drives 


through its regular patrol route. After being released from 





the car near an area of an illegal activity, the BP agent 








tries to track illegal entrants and capture them. However, 





due to the short sensor range, it usually loses contact or 








never gains contact with enemy. Then the BP releases its 
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child squad, the MiniUAV, to re-establish contact with 
enemy. With the help of the MiniUAV, the BP tracks the 


illegal entrants and finally captures them. 





There are two BP agents in each BP Car, and the BP is 





the only agent that has the ability of capturing illegal 








entrants. The BP can get information directly from the 





MiniUAV or indirectly from other sensors via the 


Communication Center. 


8. Scouts 





There are seven scouts distributed along the border to 


provide information about BP activities to the illegal 





entrants. They are stationary in their observation posts 


and invisible to the other sensors. They also extend their 





sensor range by using binoculars. Scouts communicate with 





smugglers and illegal immigrants whenever there is a 


detectable movement along the border. 











9. Illegal Immigrants 
There are 20 illegal immigrants trying to cross the 
border in the base model. Initially, they all spread out 


along the Mexican part of the border at random. Illegal 


immigrants have a treat level of two which makes them a 





less valuable target for BP and other agents. They follow 
their waypoints into U.S. to the North and try to avoid the 


BP along their way. Immigrants have no communication 








capability with other agents but they can get information 








on the BP activities via Scouts. Once illegal immigrants 








are captured by a BP agent they are disregarded by the 


model. 
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10. Smugglers 


There are two smugglers trying to cross the border by 


using illegal immigrants as distractions to the BP in the 





base model. They are randomly placed along the border and 


are given the same routes as the immigrants. Having a 











threat level of three, they are the most valuable targets 





for BP and other friendly agents. Smugglers have the 


capability of communicating with Scouts for the BP actions 











along the border. As with immigrants, they try to cross the 


border by avoiding the BP agents. 


11. Neutrals 


Neutrals can be detected by both enemy and friend and 





cause distraction until they are classified as neutral. 


They are distributed around the battlefield and given some 





waypoints to ensure they cross the border back and forth 
during the simulation. Neutrals are the agents that are 


used to create false alarms for the sensors. 
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IV. ANALYSIS 


A. MEASURES OF EFFECTIVENES 


The ultimate goal of the border protection system, 
with all of its sensors and manpower, is to deter and 


prevent illegal entries through the border. However, only 





using state of the art sensors that can detect and classify 
illegal activities are not enough to ensure protection. 


These sensors can detect and classify an activity, but 





without a tangible action, there is nothing to prevent 
illegal entrants. In the current border security system, 


this tangible action is provided by BP agents. Therefore, 











th percentag of captured illegal entrants is a 


considerable measure of effectiveness (MOE) of the overall 








system. 


There are two types of illegal entrants in the model, 


illegal immigrants and smugglers. Smugglers are high-valued 





targets EOr BP agents. Therefore, to measure the 





ffectiveness of the system, it is important to explore how 





well the system does on distinguishing between two types of 








illegal entrants. In other words, can the system focus on 





the targets of the highest value to prevent ineffective us 


of resources? 





While the percentage of captured illegal entrants is a 


good measure of overall system, due to the fact that BP 





agents are the only entities that have the ability to 
capture, this MOE is highly dependent on the number and 





capabilities of the BP agents. Therefore, the number of 
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illegal entrants classified is a more direct MOE to 





investigate th ffectiveness of the sensors. 


In the light of these ideas, the MOEs that are used in 


this study are as follows: 


e Total number of illegal immigrants captured 





e Total number of smugglers captured 


e Total number of illegal immigrants classified 


























e Total number of smugglers classified 


B. DESIGN OF EXPERIMENTS (DOE) 
i Overview 


DOE plays an integral role in the conduct of the 
Simulation study. The simulation study may use both 
gualitative and quantitative factors. Each factor can be 
set to two or more factor levels. A design point consists 


of the specification of each factor in the study. Also, 





different pseudo-random numbers (PRNs) are used to simulate 








th sam scenario in stochastic simulations. Each run of 








the simulation scenario with a specific PRN is called a 


replicat (Kleijnenet al., 2005). 





2. Important Factors and Ranges 


There are numerous parameters in the MANA model that 





can be varied to explore the effects on the outcome. Due to 


computational limitations, a subset of these parameters, 





that is deemed to be most likely to cause major effects on 








th selected MOES, is chosen as factors for the 
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experimental design. 


then 


determined. 





extended to determin 


Table 10 shows 


Current 


th 


capabilities 











in the experimental design. 


of the 





Categorical Factor 


With MiniUAV 





Without 
MiniUAV 





agents 


Reasonable ranges of these factors are 


are 


maximum ranges for the factors. 


the factors and their ranges that are used 





Ranges of Numeric Factors 





























































































































Low Low High High 
Level Level Level Level 
(MANA (Real (MANA (Real 
Factor Name Units) World) Units) World) decimals 
Fuel (Represent 160 
1 endurance) 800 40 min. 3200 min 0 
2 Turret Height 150 m. 150 m. 400 m. | 400 m. 0 
Probability of 
3 | MiniUAV | Classification (max) 0.3 0.3 1 1 2 
Classification 1000 
4 Range (max) 4 200 m 20 m. 0 
1500 
5 Detection Range 6 300 m 30 m. 0 
60 
6 Comm. Latency 0 O sec 20 sec 0 
MinivUAV 
7 Reliability 60 60% 100 100% 0 
8 Turret Height 20 20 m. 50 50 m 0 
Probability of 
9 Classification (max) 0.3 0.3 1 al 2 
Tower Classification 10000 
10 Range (max) 50 2500 m. 200 m. 0 
15000 
a Detection Range 75 3750 m. 300 m. 0 
Probability of 
12 Classification (max) 0.3 0:23 il 1 2 
Predator | Classification 10000 
13 B Range (max) 50 2500 m. 200 m. 0 
15000 
14 Detection Range 75 3750 m. 300 m. 0 
60 
TS: Comm. Latency 0 O sec. 20 sec 0 
Center 
16 Reliability 60 60% 100 100% 0 
BP Car 
17 | (Package) of BP car agents 2 2 5 5 0 
Table 10. Important Factors and Ranges 
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3. NOLH Design 


Minimum and maximum levels of these factors form the 





frame of the experimental region. To explore the main 
effects and the interactions of the factors, we need to 


select points to sample within that experimental region. 











Despite the technological improvements on computers, it is 
still practically impossible to run a model for all 


possible points in an experimental region. For example, if 





we have 10 factors with 10 levels each, we need 10'° = 10 








billion design points for a factorial design. Therefore, we 
need sophisticated techniques to find an efficient number 
of design points, which together allow for maximum 


information to be gained from the experiment. 


According to Sanchez, Latin hypercube (LH) designs 
provide a flexible way of constructing efficient designs 
for quantitative factors and have some of the space-filling 
properties of factorial designs, but require orders of 


magnitude less sampling (Sanchez, 2008). 


Figure 14 shows the excellent space filling behavior 





of the Nearly Orthogonal Latin Hypercube (NOLH) design, 





with less sampling requirement, compared to factorial 


designs. 




















24 factorial: 16 design points 44 factorial: 256 design points NOLH: 17 design points NOLH: 257 design points 


Figure 14. Scatterplot Matrices for Selected Factorial and NOLH 
Designs. (From: Sanchez, 2008) 
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An Excel spreadsheet, which was developed by Sanchez, 





is used to construct the NOLH experimental design for this 
study (Sanchez, NOLHdesigns spreadsheet, 2005). This 
spreadsheet provides five different NOLH designs’ for 
various numbers of factors. The first four designs, for up 


to 22 factors, are from Tom Cioppa's 2002 PhD dissertation 





(2002). The last one, for experiments with up to 29 


factors, is provided by Andy Hernandez. 


According to Table 10, 17 factors are chosen for the 





experimental design. However, classification ranges are set 
to 2/3 of the detection ranges in the model which makes 


these pairs linear combinations of each other. Therefore, 





the number of independent factors used in the experimental 





design was 14. Minimum and maximum values of these factors 
were then typed into the 29-factor spreadsheet in Sanchez’s 


NOLH design spreadsheet. As a result, 257 levels of these 





14 factors are acquired automatically. The reason for using 





the 29-factor design spreadsheet, despit having only 14 





factors, is th presenc of categorical factor, “with 


MiniUAV?.” The inclusion of this binary factor, met the 





desire to have more design points, in order to reduce the 


maximum pairwise correlation and to increase the degrees of 





freedom for error. 


The first seven rows of the NOLH design, constructed 
for this study, are shown in Figure 15 (See Appendix A for 


the whole design). There are additional columns for 





dependent factors in this experimental design spreadsheet. 








Three of them are the detection ranges for the sensors as 


mentioned above. 


5:3 


4 0 60 
20 20 100 


| MiniUAVTurretHt 
o 2 = MiniUAVClassRng 5 
MiniVAVCommLat 5 
J2 23838 BS MiniVAvcommRel 5 
‘wrClassRng(2) 
‘wrP Class(1) 
PredClassRng(1) 
3 u CommCentLat 


co  S MiniVAVDetRng 
WithUAV?(Y/N) 


Factor Name 
MiniUAVFuel 
TwrClassRng 
TwrDetRng 
TwrClassRng(1) 
PredDetRng(2) 
& PredClassRng(2) 
BPCar5(Sqd9) 
BPCar1(Sqd7) 
BPCar4(Sqd8) 


TwrPClass 


‘© TwrPClass(2) 
& PredClassRng 
& PredDetRng 

& PredPClass(2) 
~ PredPClass(1) 


o 
i PredPClass 


o 


0.88 177 266 
0.86 184 276 
0.82 192 288 
0.73 153 230 0. 102 0.95 
0.71 193 290 0. 129 0.94 
0.74 138 207 0.89 138 92 0.95 


223 
283 
38 
Bas 
3 
i 
— 
z 
R 
np 


0No No No No 
1 Yes Yes Yes Yes 
0 Yes No No No 
1 Yes Yes Yes Yes 
0 Yes No_No_No 


o 
& 
Rg 
R 
BQ x S PredDetRng(1) 


= 


0.31 132 198 
0.25 172 258 
0.71 119 179 
_0.67_ 138 207 


1775 190 
1681 290 
1859 242 
903 303 





NNO 
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Figure 15. Experimental Design Spreadsheet 


There are also two sets of additional columns for 








sensor ranges and corresponding probability of 





classification values. While modeling Towers and Predator B 
in MANA, “advanced mode” is used for their sensors. What 


this means is that instead of giving only one 





classification range and thereby creating a 
“cookie-cutter” sensor, three levels of ranges and 
corresponding probabilities are used in decreasing order. 


Given a maximum classification range, medium and minimum 





ranges are calculated by multiplying max range by 2/3 and 


1/3. Similarly, values for probability of classification 





are calculated for medium and minimum ranges by using 
P(Class) of the max range as an input. Formulas, used for 
medium and minimum ranges, are 1l-exp(-3*P(Class for max 


range)) and 1l-exp(-4*P(Class for max range)). 


Finally, the last four columns describe the binary 
“WithUAV” variable as well as the specification of which 
specific BP cars are used in the model, corresponding to 


the various levels of the “Number of BP Cars” variable. 


To see the space filling property of the NOLH design, 


a scatter plot is presented in Figure 16 that shows the 
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pairwise projections of the levels of the factors used in 


the experimental design. 
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Figure 16. Scatterplot Matrix for the Experimental Design 
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Cc. MODEL RUN 
1. Limitations of MANA 


Although MANA has some simple tools for analysis, it 





is not capable of constructing a large scale experimental 





design and does not provide sufficient analysis capability. 


Therefore, Sanchez’s NOLH design spreadsheet is used to 





construct th xperimental design. Another piece of 
software, called XStudy, is used to create MANA scenario 
files, one corresponding to each design point. For every 
design point, 30 replications are run using the random seed 


generator of MANA, in order to see the variability within 





each scenario. 


2. XStudy 


The XStudy software is used to specify the number of 
random replications and the mapping of each column in the 
design to a particular element in the MANA scenario file. 


As a result, a study design file, in xml format, is 





generated. Next, software called OldMcData generates th 
individual MANA scenario files for each design point, using 
the design spreadsheet, the base case MANA scenario xml 


file, and the xstudy xml file. As the last step, software 





called Condor manages the distribution and collection of 
individual jobs across a set of processors. Condor is open 
source software available from the University of Wisconsin 


at http://www.cs.wisc.edu/condor. After the jobs have 








completed, OldMcData contains a postprocessor to gather all 





model outputs into one file, along with the design point 


values. XStudy and OldMcData were developed by SEED Center 
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Research Associate Stephen Upton, and are available for 


download at the SEED Center web site at 





http://harvest.nps.edu. 





D. STATISTICAL TOOLS AND ANALYSIS 
1. JMP 


After running the model for 257 design points and 30 


replications of each, 7710 data points are acquired for 





each four MOEs. JMP is the primary tool used to analyze the 


results. 


According to the JUMP website, JUMP is a dynamic, 


interactive, visual, and easy to use statistical software 





package developed by SAS Institute, Inc. In JMP, commands 








are available to interactively extend analyses and further 
explore the results. The rows in JUMP data tables are 
dynamically linked which allows the user to identify any 


data point in graphs easily. 


JMP has a Custom designer tool to construct 


experimental design, supporting screening, Response Surface 








Methodology (RSM), and mixtur xperiments. JUMP also has 
the main fitting tools such as linear fitting, mean- 
dispersion models, partition trees, nonlinear regression, 
and Neural Nets, supported by visual interfaces 


(Creightonet al., 2005). 


2. Statistical Methods Used To Analyze Results 


In order to fit models, the 7710 data points are 
summarized by their means and standard deviations, thus 
collapsing the data set for analysis down to 257 data 


points. The distribution of these means for each MOE is 
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explored graphically and with basic summary statistics. The 


data are also examined for outliers or any type of anomaly. 





To further explore the effects of the categorical 
factor, MiniUAV, a comparison test is conducted for each 


MOE. The T-Test is used to reveal any significant 





differenc between the means of MOEFS with MiniUAV and 
without MiniUAV. 





Linear stepwis regression models are fitted to 





explore the effects of various factors on different MOEs. 





For each MOE, the regression models are fit separately for 


the subset of points where the MiniUAVS are in the 





scenario, and then fit for the subset of points where the 


MiniUAVs are not in the scenario. The difference between 





the two is that the UAV performance factors are considered 


for inclusion in the model only when the MiniUAVs are in 





the scenario. Main effects and two-way interaction terms 





of the factors are considered for inclusion in these 





models. 


Finally, regression trees are used to discover the 


threshold values of the factors that are found to be 





important in linear regression models. 
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Vv. RESULTS 


As SUMMARY OF STATISTICS 
alk Total Number of Captured Immigrants 


There are 20 immigrants that are trying to cross the 
border during the simulation. They can only be captured by 
BP agents. Figure 17 shows the distribution of the number 
of captured immigrants and the summary statistics. Although 


the number of BP agents is varied between four and ten, the 





mean of the number of captured immigrants is 3.04 for the 
“without MiniUAV” scenarios and 5.49 for the “with MiniUAV” 
scenarios. Despite the slight improvement in mean by using 
the MiniUAV, these mean values are still low compared with 


the number of immigrants. 





Total Number of Captured Immigrants | Total Number of Captured Immigrants 
without MiniUAV(Mean) with MiniUAV(Mean) 





= 
fo) 












no oOo fF aA DBD NN O © 


1 
































Moments Moments 

Mean 3.0304688 Mean 5.4937984 
Std Dev 1.0914607 Std Dev 1.89234 
Std Err Mean 0.0964724 Std Err Mean 0.1666113 
upper 95% Mean 3.2213702 upper 95% Mean 5.8234673 
lower 95% Mean 2.8395673 lower 95% Mean 5.1641296 
N 128 N 129 

Figure 17. Distribution of Total Number of Captured Immigrants 
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2. Total Number of Captured Smugglers 





There are only two smugglers within 20 illegal 
immigrants that are trying to cross the border during the 


Simulation. They are high priority targets for the sensors 





and BP agents, and they can only be captured by BP agents. 
Figure 18 shows the distribution of the number of captured 
immigrants and the summary statistics. The distribution of 
the captured smugglers has some gaps, due to the small 
number of smugglers used. The mean of the total number of 
captured smugglers is 0.3 for the “without MiniUAV” 
scenarios and 0.51 for the “with MiniUAV” scenarios. 


Similar to the number of captured immigrants, the 





effectiveness of BP agents in capturing smugglers is not 


satisfactory in both scenarios. 





Total Number of Captured Smugglers | Total Number of Captured Smugglers 
without MiniUAV(Mean) with MiniUAV(Mean) 




















Moments Moments 

Mean 0.2986979 Mean 0.5129199 
Std Dev 0.1023266 Std Dev 0.1566971 
Std Err Mean 0.0090445 Std Err Mean 0.0137964 
upper 95% Mean 0.3165953 upper 95% Mean 0.5402185 
lower 95% Mean 0.2808005 lower 95% Mean 0.4856213 
N 128 N 129 

Figure 18. Distribution of Total Number of Captured Smugglers 


60 


3. 


Because 


illegal entrants, 


only BP 


agents 


are 


capable 


Total Number of Classified Immigrants 


of 


the total number of classified 


capturing 


immigrants 





is examined to see the effectiveness of the sensors more 


independently. 
of 
and 19.02 for 


number 


MiniUAV” 


result means 


immigrant 


once. 





classified 


immigrants 





that, 





in both 


“with MiniUAV” 


cases, 


With a slight difference, 


is 


mean of 
18.9 for 
(Figure 


almost 





the total 
“without 
19). This 


every illegal 





is classified with one of the sensors at least 





Mean(TotimmigClassified) 











Mean(TotimmigClassified) 














4. 


Similar to the case with 
scenarios all of the smugglers are cl 


the sensors at least once. 





Immigrants 


illegal 


Total Number of Classified Smugglers 


immigrants, 
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assified wi 


Moments Moments 
Mean 19.020672 
M 18.905469 
Sid Dev 0.9787039 Std Dev 0.9481634 
Std Err Mean 0.086506 Std Err Mean 0.083481 1 
upper 95% Mean 19.076649 upper 95% Mean 19.185854 
lower 95% Mean 18.734289 lower 95% Mean 18.85549 
N 128 N 129 
Figure 19. Distribution of Total Number of Classified 


in most 


th one of 





Figure 20 shows that 





mean of 


the total number of classified smugglers is roughly 1.93 


whether the MiniUAV is used or not. 





Mean(TotSmugClassified) Mean(TotSmugClassified) 


















































Moments 














Moments Mean 1.9341085 
Mean 1.9338542 Std Dev 0.0893131 
Std Dev 0.0814334 Std Err Mean 0.0078636 
Std Err Mean 0.0071978 upper 95% Mean 1.949668 
Upper ee a Mean ee lower 95% Mean 1.9185491 
lower 95% Mean 1.9196111 N 429 
N 128 
Figure 20. Distribution of Total Number of Classified 





Immigrants 


B. T-TEST FOR COMPARISON OF SAMPLE MEANS 


To answer the second research question, which is how 
the UAV can contribute to border protection, scenarios 


without MiniUAV” and “with MiniUAV” are compared in terms 





of their corresponding MOEs. 


1. Total Number of Captured Immigrants 


According to the t-Test for the means of total number 


of captured immigrants, there is a significant difference 








in 95% confidenc level between the scenarios “without 





MiniUAV” and “with MiniUAV”. It is seen in Figure 21 that 
using the MiniUAV has a positive effect on the mean of the 


total number of captured immigrants. 
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Oneway Analysis of Mean of Captured t Test 


Immigrants Yes-No 
Assuming unequal variances 








95 Difference 2.46333 t Ratio 12.79479 

3 Std Err Dif 0.19253 DF 204.9936 
a Upper CL Dif 2.84291 Prob > |t| <.0001 
g? Lower CL Dif 2.08374 Prob >t <.0001 
8 6 Confidence 0.95 Prob <t 1.0000 
ne} 
BS 
S 4 
= 

34 

274 











WithUAV 























Figure 21. Comparison Test for Total Number of Captured 
Immigrants 


2. Total Number of Captured Smugglers 


The comparison test yields the same overall result for 


the total number of captured smugglers. We are 95% 





confident that there is a significant difference between 
the scenarios “without MiniUAV” and “with MiniUAV” (See 


Figure 22). 











Oneway Analysis of Mean of Captured t Test 
Smugglers Yes-No 
14 Assuming unequal variances 
a 
a P Difference 0.214222 t Ratio 12.98569 
0.9 Std Err Dif 0.016497 DF 220.5971 
@ 0.8 i Upper CL Dif 0.246733 Prob > <.0001 
O 0.7 H It 
N o6 H Lower CL Dif 0.181711 Prob >t <.0001 
O05 : H Confidence 0.95 Prob <t 1.0000 








No Yes 
WithUAV 


























Figure 22. Comparison Test for Total Number of Captured 
Smugglers 
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3. Total Number of Classified Immigrants 


Although there is significant evidence of the 
contribution of the MiniUAV in capturing illegal 


immigrants, there is no significant difference between the 








means of the total number classified immigrants. For most 


scenarios, the sensors classify all of the immigrants, and 





the distribution of the total number of classified 
immigrants is highly skewed towards the upper end. Thus 
there is very little variability between the two scenarios. 


Increasing the number of illegal immigrants or their 





concealment parameter would make it hard for sensors to 


classify the immigrants, and may produce different results. 











Oneway Analysis of Mean(TotimmigClassified) | t Test 
By WithUAV Yes-No 
; hd Ps a ici _ “ 
B19 | , 
= Difference 0.11520 t Ratio 0.958284 
A Std Err Dif 0.12022 DF 254.6028 
ra 18 Upper CL Dif 0.35195 Prob > |t| 0.3388 
aD Lower CL Dif -0.12155 Prob >t 0.1694 
E 17 Confidence 0.95 Prob <t 0.8306 
6 
Ee 
5 16 
®o 
= 
15 
aa nee 0.4 -0.3 0.1 0.0 0.1 0.2 0.3 0 
WithUAV aad earns 
























Figure 23. Comparison Test for Total Number of Classified 
Immigrants 





Due to the lack of information gained from examining 





the total number of classified immigrants, another measure 





of effectiveness, the time to classify, was sought. Thus, 
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if all of the immigrants are classified through the 


Simulation, how long does it take? 


The mean time to classify all of the 





immigrants was 


7028 time steps, 5.86 hours, for the “without MiniUAV” and 


6846 time steps, 5.7 hours, for the “with MiniUAV” scenario 


(See Figure 24). According to the comparison 
a significant decrease in mean time to 


immigrants when the MiniUAV is used. 


test, there is 


classify all 





Oneway Analysis of MaxTimelmmigClassified By | Missing Rows 
WithUAV one 

t Test 
Yes-No 





Difference 
Std Err Dif 
Upper CL Dif 
Lower CL Dif 
Confidence 


Assuming unequal variances 


-181.26 t Ratio -3.54748 
51.10 DF 3885.238 
-81.08 Prob > |t| 0.0004 

-281.44 Prob >t 0.9998 

0.95 Prob <t 0.0002 





MaxTimelmmigClassified 








No Yes 
WithUAV 











Momenis(without MiniUAV) 


Mean 7028.1598 








-200 -100 -50 0 50 100 150 20 





Momenis(with MiniUAV) 





Mean 6846.898 

moa 1585.432 Std Dev 1601.7441 

Std Err Mean 36.3531 76 Std Err Mean 35.905965 

upper 95% Mean 7099.4561 upper 95% Mean 6917.3152 

lower 95% Mean 6956.8635 lower 95% Mean 6776.4807 
Figure 24. Comparison Test for the Time All Immigrants are 





Captured 
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4. Total Number of Classified Smugglers 


The result of the comparison test for the total number 
of classified smugglers is very similar to immigrants. 


There is no significant difference between the means of 





classified smugglers in terms of MiniUAV usage (Figure 25). 
Due to the similarly skewed distribution of classified 
smugglers, the same conclusion can be made as with the 


immigrants. 





Oneway Analysis of Mean(TotSmugClassified) t Test 








By WithUAV Yes-No 

Ss 
2 1.9 Difference 0.00025 t Ratio 0.02386 
2 Std Err Dif 0.01066 DF 253.1998 
ra Upper CL Dif 0.02125 Prob > 0.9810 
Dm |t| 
28 Lower CLDif 0.02074 Prob >t 0.4905 
g Confidence 0.95 Prob <t 0.5095 
kr 
se 1.7 
o 
® 
= 

1.6 

No Yes 
WithUAV 


























Figure 25. Comparison Test for Total Number of Classified 
Smugglers 


However, there is a significant difference between the 





means of the time to classify all of the smugglers. The 





mean time to classify all of the smugglers is 4398 time 





steps, 3.67 hours, for the “without MiniUAV” and 4228 time 





steps, 3.52 hours, for the “with MiniUAV” scenario (Figure 


26). 
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Oneway Analysis of MaxTimeSmugClassified By Missing Rows 
WithUAV Bas 
11000 tTest 
Yes-No 
10000 
8 9000 Assuming unequal variances 
‘® 8000 Difference -169.70 t Ratio -4.60733 
8 7000 Std Err Dif 36.83 DF 7211.35 
2 Upper CL Dif -97.50 Prob > <.0001 
= 6000 | 
® 5000 Lower CL Dif -241.90 Prob >t 1.0000 
£ Confidence 0.95 Prob <t <.0001 
FE 4000 
& 3000 
2000 
1000 
No Yes 
WithUAV 
-200 -100 -50 O 50 100 150 20 
Moments(with MiniUAV : saa 
Oments(wimnout ve Moments(with MiniUAV) 
Mean 4398.3201 
Std Dev 1589.4047 Mean 4228.6234 
Std Err Mean 26.493758 Std Dev 1540.5122 
upper 95% Mean 4450.2644 Std Err Mean 25.586514 
lower 95% Mean 4346.3758 upper 95% Mean 4278.7888 
N 3599 lower 95% Mean 4178.458 
N 3625 
Figure 26. Comparison Test for the Time All Smugglers are 
Captured 
CG. LINEAR REGRESSION 








A linear multiple stepwise regression model is fitted 





to examine the ffects of the factors on the MOEs. As 


stated previously, these models are fit in two stages for 





each MOR - first selecting all of the points “with MiniUAV” 





(and considering the UAV performance factors for 


inclusion), and then selecting all of the points “without 








MiniUAV”. All factors that are used in the experimental 





design and their two-way interaction terms are considered 


for inclusion in the models. 
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1. Total Number of Captured Immigrants 
a. Without MiniUAV 


The p-value of the linear regression model for 
the number of captured immigrants is less than 0.0001, 
meaning that there is a Significant linear relation between 


the response and at least one of the regressors. The 








RSquare value is 0.91, which indicates that the regressors 
explain successfully the variability of the response, in 


this case the total number of captured immigrants. 


There are six important terms that affect the 





response at the 0.95 confidence level (Table 11). The most 
important one is the “Number of BP Cars.” This is very 
meaningful, because BP agents are the only agents that are 


capable of capturing illegal immigrants. Because there ar 








two BP agents in each BP car, the number of BP cars 








directly affects the number of BP agents. 


Another highly important factor is “Communication 
Center Reliability”. A perfect (100%) reliability for 
communication systems means that a communication message 
(containing contact location information) from one agent to 


another is always successfully delivered. Becaus th 








Communication Center is the most important element of the 





communication system, civil is not surprising to see 
“Communication Center Reliability” as the second important 


factor. 


“Tower Classification Range” shows up as another 








important factor. There are also some interactions that are 














important. Interactions of “Communication Center Reliance” 


with “Tower Height” and also with “Predator Probability of 





Classification” are two of them. Finally “Tower 
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Classification Range” and “Tower Probability of 














Classification” is another important interaction. 





Response Mean( Sqd23Cas) 
Actual by Predicted Plot 
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oa 


Mean( Sqd23Cas) 
Actual 





1 2 3 4 5 6 
Mean( Sqd23Cas) Predicted 
P<.0001 RSq=0.91 RMSE=0.3434 








Summary of Fit 


RSquare 0.910388 

RSquare Adj 0.901037 

Root Mean Square Error 0.343356 

Mean of Response 3.030469 

Observations (or Sum Wgts) 128 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 12 137.73569 11.4780 97.3592 
Error 115 13.55771 0.1179 Prob > F 
C. Total 127 151.29339 <.0001 


Sorted Parameter Estimates 
Term 
NumBPCars 
CommCentRel 
TwrClassRng 
(TwrHt-35.1875)*(CommCentRel-63.5547) 
(TwrPClass-0.54305)*(TwrClassRng-126.367) 
(PredPClass-0.54633)*(CommCentRel-63.5547) 















Estimate | Std Error|t Ratio|t Ratio 


|_0.969351| 0.041503] 23.36 | {Tn 
[0.0086565[ 0.001515] 5.71 |[{ 
| -0.003402[ 0.001022] -3.33|{ 
| -0.00068[ 0.00022] -3.08|[ TT 
[ -0.011695] 0.004585] -2.55|[ 

00 | 


0.0137206 | 0.006856 





Prob>|t| 




























Table 11. Linear Regression Results for Total Number of 
Captured Immigrants without MiniUAV 
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b. With MiniUAV 


With a p-value of less than 0.0001 and RSquare 


value of 0.965, the regressors explain successfully the 





variability of the response for the scenarios with the 


MiniUAV. 


Now that the factors related to the performance 
of the MiniUAV are added to the model there are some 
differences in the important factors. “Number of BP Cars” 
is still the most important factor, though with a slight 


difference from the “without MiniUAV” scenario. This factor 





is more important in this model because increasing the 





number of BP cars also increases the number of MiniUAV 


used. 


Two new important factors related with MiniUAV 


are “MiniUAV Classification Range” and “MiniUAV Fuel”, 





which designates the endurance of the MiniUAV. 





Interaction between probability of classification 





of Tower and Predator is the fourth important factor. See 
Table 12 for the other less important factors and 


interactions. 
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Response Mean( Sqd23Cas) 
Actual by Predicted Plot 
10 





Mean( Sqd23Cas) 
Actual 





- NWO FPUAHDN WO 


12 3 4 5 6 7 8 9 10 
Mean( Sqd23Cas) Predicted 
P<.0001 RSq=0.96 RMSE=0.4199 








Summary of Fit 


RSquare 0.965 

RSquare Adj 0.950769 

Root Mean Square Error 0.419875 

Mean of Response 5.493798 

Observations (or Sum Wgts) 129 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 37 442.31890 11.9546 67.8102 
Error 91 16.04281 0.1763 Prob > F 
C. Total 128 458.36171 <.0001 


Sorted Parameter Estimates 





Term Estimate Std t/t Ratio Prob>|t| 
Error| Ratio 

NumBPCars |1.6287062| 0.055248| 29.48|[ | <.0001 

MiniUAVClassRng [0.1122961 | 0.008246| 13.62/[[ Te id <.0001 

MiniUAVFuel 10.0007129| 6.734e-5] 10.59/— [jm _]| <.0001 


(TwrPClass-0.55698)*(PredPClass-0.55372) 5.9925185| 1.09983] 5.45/- 1m | <.0001 
(MiniUAVTurretHt-275.651)*(TwrHt-34.8217) -0.000321| 9.368e-5| -3.43/— sm] 0.0009 


(PredClassRng-124.209)*(CommCentRel- -0.000139] 4.982e-5 ——— 0.0066 
66.4419) 


(MiniUAVFuel-2016.79)*(CommCentRel-66.4419) | 1.2168e-5| 4.396e-6| _2.77|- 1 | 0.0068 
(MiniUAVFuel-2016.79)*(TwrPClass-0.55698) __|0.0007813| 0.000285| 2.74/- 1m | 0.0074 


(MiniUAVCommRel-79.3256)*(NumBPCars- 0.0097734] 0.003841 | 244 0.0126 
3.51938) 


MiniUAVFuel-2016.79)*(TwrClassRng-123.651) | 4.4024e-6| 1.838e-6| 240/719 | 0.0186 








MiniUAVCommRel 10.0083524| 0.003525] 2.37/> 168 _]|_ 0.0199 
(MiniUAVFuel-2016.79)*(NumBPCars-3.51938) _[0.0001571| 6.683e-5| 2.35|/[ Td (0.0209 
(TwrHt-34.8217)*(TwrPClass-0.55698) 10.0648917| 0.029069| 2.23/— m8 _]|_ 0.0280 
TwrHt | 0.009718] 0.004499] 2.16/7 Tr —idCO0.0334 
(TwrPClass-0.55698)*(TwrClassRng-123.651) _[0.0125866| 0.006103] _2.06/- 19 ___]| 0.0420 
PredPClass | -0.385112| 0.190985] -2.02/[7 Td CO0.0467 
(PredPClass-0.55372)*(CommCentLat-10.2403) |0.0635376|0.031699| 2.00/- 15 ]| 0.0480 
Table 12. Linear Regression Results for Total Number of 


Captured Immigrants with MiniUAV 
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2. Total Number of Captured Smugglers 
a. Without MiniUAV 


Although the p-value is less than 0.0001, which 





means at least one regressor has a linear relation with the 
number of captured smugglers; the RSquare value is rather 


low, 0.36, meaning that the regressors are not enough to 





explain the variability in the response. The reason for 





that may be that there is a non-linear relation with some 
of the independent variables, or there are some other 
uncontrollable factors such as the random patterns of the 
smugglers. The only important factor that has a linear 


relation with this MOE is again the number of BP cars. 
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Response Mean( Sqd24Cas) 
Actual by Predicted Plot 


Bos al EZ 





Mean( Sqd24Cas) 


0.1 0.2 0.3 0.4 0.5 


Mean( Sqd24Cas) Predicted 
P<.0001 RSq=0.36 RMSE=0.0843 








Summary of Fit 


RSquare 0.361828 

RSquare Adj 0.318925 

Root Mean Square Error 0.084338 

Mean of Response 0.298438 

Observations (or Sum Wgts) 128 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 8 0.4799115 0.059989 8.4337 
Error 119 0.8464426 0.007113 Prob > F 
C. Total 127 1.3263542 <.0001 


Sorted Parameter Estimates 





Term Estimate | Std Error|t Ratio|t Ratio Prob>|t| 
NumBPCars 0.0566461/0.007861| 7.21/—- —_7 _| | <.0001 


(PredPClass-0.54633)*(NumBPCars-3.48438) _[0.0786461| 0.044517|_1.77/[- | mm __]| 0.0799 
(PredPClass-0.54633)*(CommGentLat-9.77344) | -0.01051| 0.006007] -1.75[—- mma 1 ___]| 0.0828 











TwrClassRng | -0.000249] 0.000172 -1.45[[ Tm Tid (0.1506 
TwrPClass }0.0531324] 0.036804] 1.44 |[ | 0.1515 
(CommCentLat-9.77344)*(NumBPCars-3.48438) | -0.002017| 0.001444|_-1.40|[ TT | (0.1650 
PredPClass [0.0409693| 0.036936[ 1.11/71 a 1 ___]|_ 0.2696 
CommCentLat 0.0013002| 0.001281] 1.01/- | m1 | 0.3122 
Table 13. Linear Regression Results for Total Number of 


Captured Smugglers without MiniUAV 


b. With MiniUAV 


The variability in the number of § captured 


smugglers is better explained in “with MiniUAV” scenarios, 





with an RSqaure value of 0.65. The important factors other 





than the number of BP cars are “MiniUAV Classification 





Range”, “MiniUAV Fuel” and interaction between probability 
13 





Table 14. 


of classification of the Tower and Predator, which is very 





Similar to the model for the number of captured immigrants. 


See Table 14 for the 
interactions. 


Response Mean( Sqd24Cas) 
Actual by Predicted Plot 


other less important factors and 





Mean( Sqd24Cas) 


Mean( Sqd24Cas) Predicted 





P<.0001 RSq=0.65 RMSE=0.1013 





, 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 





Summary of Fit 





RSquare 0.654576 

RSquare Adj 0.582884 

Root Mean Square Error 0.101296 

Mean of Response 0.512145 

Observations (or Sum Wgts) 129 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 

Model 22 2.0610968 0.093686 9.1304 

Error 106 1.0876543 0.010261 Prob > F 

C. Total 128 3.1487511 <.0001 

Sorted Parameter Estimates 

Term Estimate Std t|t Ratio Prob>|t| 
Error| Ratio 

NumBPCars 10.0669452| 0.009713 6.89/- |_| | <.0001 

MiniUAVClassRng [0.0110488] 0.001968] 5.61|[ | dn | <.0001 

MiniUAVFuel [5.8111e-5| 1.325e-5| 4.39/[ ed <.0001 


(MiniUAVClassRng-11.8682)*(PredPClass- 
0.55372) 
(TwrHt-34.8217)*(PredPClass-0.55372) 
(MiniUAVClassRng-11.8682)*(TwrHt-34.8217) 
(MiniUAVFuel-2016.79)*(MiniUAVClassRng- 
11.8682) 

CommCentRel 
(MiniUAVClassRng-11.8682)*(CommCentRel- 
66.4419) 
(PredPClass-0.55372)*(CommCentRel-66.4419) 
(MiniUAVFuel-2016.79)*(PredPClass-0.55372) 
(MiniUAVClassRng-11.8682)*(PredClassRng- 
124.209) 


-0.022985| 0.007644| -3.01/- _ [1 ]| 0.0033 
0.0006875| 0.00025] 2.75|[ Tr sid, SC(0.0071 


0.0011936| 0.000466[ 2.56|[ Tn ——id|CO.0118 


-0.005364| 0.002498] -2.15/[— mmm 1 | 0.0340 
-0.000136| 6.579e-5| -2.07/— sm ~~ ~——S—S—~dY:Cé«.0408 
0.0001034| 5.055e-5|  2.04/— 1 mam | :~«(0.0433 








Linear Regression Results for Total Number of 


Captured Smugglers with MiniUAV 
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3. Total Number of Classified Immigrants 
a. Without MiniUAV 


The model for the number of classified immigrants 








can explain most of the variability with an RSquare value 





of 0.93. Because this MOE measures th ffectiveness of the 





sensors, it is to be expected not to see the number of BP 





cars aS an important factor in this model. 


The two most important factors for this model are 


“Predator Classification Range”, “Tower Classification 











Range” and the interaction term of these two. See Table 15 


for the other important terms and interactions. 


75 


Actual by Predicted Plot 





Mean(TotlmmigClassified) 
Actual 





15 16 17 18 19 20 


Mean(TotlmmigClassified) Predicted 
P<.0001 RSq=0.93 RMSE=0.2833 








Summary of Fit 


RSquare 0.925451 

RSquare Adj 0.916214 

Root Mean Square Error 0.283293 

Mean of Response 18.90547 

Observations (or Sum Wgts) 128 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 14 112.57959 8.04140 100.1982 
Error 113 9.06881 0.08025 Prob > F 
C. Total 127 121.64839 <.0001 


Sorted Parameter Estimates 











Term Estimate | Std Error|t Ratio|t Ratio Prob>|t| 
PredClassRng [0.0149257| 0.000589] 25.35|(" | [im |<.0001 
TwrClassRng 10.0124413] 0.000583| 21.34[— 1 jms | <.0001 
(TwrClassRng-126.367)*(PredClassRng-125.805) | -0.000185| 1.467e-5| -12.63|[ mmm —~d|C<.0001 
TwrHt 10.0147774| 0.002913] 5.07/[— 1mm —__]| <.0001 
(TwrHt-35.1875)*(PredClassRng-125.805) [ -0.000319| 7.431e-5| -4.30[[ mm id <.0001 
TwrPClass 10.3688509| 0.123735] 2.98/[—- 1m | 0.0035 
PredPClass 10.3487601| 0.123999] 2.81|[ | (0.0058 
(TwrPClass-0.54305)*(PredPClass-0.54633) 1.5661243| 0.689581|  2.27|/- 1m] 0.0250 
Table 15. Linear Regression Results for Total Number of 


Classified Immigrants without MiniUAV 





b. With MiniUAV 

The model for the total number of classified 
immigrants with MiniUAV is very similar to the model 
without MiniUAV with an RSquare value of 0.97. The first 
three important factors are exactly the same, and there are 
some other factors and interactions added to the model, 


related to the performance of the MiniUAV (See Table 16). 
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Response Mean(TotimmigClassified) 
Actual by Predicted Plot 
20 





19 


18 


17 


16 


Mean(TotlmmigClassified) 
Actual 


15 





15 16 17 18 19 20 


Mean(TotlmmigClassified) Predicted 
P<.0001 RSq=0.97 RMSE=0.2183 








Summary of Fit 


RSquare 0.966884 

RSquare Adj 0.947014 

Root Mean Square Error 0.218255 

Mean of Response 19.02067 

Observations (or Sum Wgts) 129 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 48 111.26293 2.31798 48.6608 
Error 80 3.81083 0.04764 Prob > F 
C. Total 128 115.07376 <.0001 


Sorted Parameter Estimates 








Term Estimate | Std Error x t Ratio 

PredClassRng 0.0125466 | 0.000646 =r 

TwrClassRng 0.011258] 0.000734 [19.34 
(TwrClassRng-123.651)*(PredClassRng- -0.000188] 1.726e-5| - 

124.209) 

(MiniUAVCommLat-10.1783)*(PredClassRng- | 0.0006386] 0.00011 

124.209) 

(MiniUAVTurretHt-275.651)*(MiniUAVCommLat- | -0.000275| 5.325e-5 ——s 
10.1783) 

TwrHt 0.0172466| 0.003501] 4.93|- Tm | 
(MiniUAVTurretHt-275.651)*(NumBPCars- 0.0020441 | 0.000451 | 
3.51938) 

CommCentRel 0.0048627| 0.001075| _4.52/- fmm | 
MiniUAVClassRng 0.019402/ 0.004467 4.34/- 1mm || 


(MiniUAVCommRel-79.3256)*(PredPClass- -0.048294| 0.011251 | 429 | 
0.55372) 
(MiniUAVClassRng-11.8682)*(CommCentLat- | 0.0041055] 0.000971 | 428 
10.2403) 
(MiniUAVCommRel-79.3256)*(TwrClassRng- 0.0002207] 5.936e-5 | 
123.651) 


MiniUAVFuel 0.0001132| 3.143e-5| 3.60/- 1mm || 


(MiniUAVCommLat-10.1783)*(NumBPCars- 0.0153324] 0.004492 | 341 
3.51938) 


(MiniUAVCommRel-79.3256)*(TwrHt-34.8217) _|0.0010605| 0.000316] 3.35|(" 1m || 


(MiniUAVTurretHt-275.651)*(CommCentLat- 0.0001851] 5.574e-5 —_—_—_— | 
10.2403) 


(PredPClass-0.55372)*(CommCentRel-66.4419) | -0.016376| 0.004949 -3.31/——_ 7 | 
(TwrPClass-0.55698)*(NumBPCars-3.51938) -0.474212| 0.144842[ -3.27/— m1] | 
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Prob>|t| 
<.0001 
<.0001 
<.0001 
<.0001 
<.0001 


<.0001 
<.0001 


<.0001 
<.0001 
<.0001 
<.0001 
0.0004 


0.0005 
0.0010 


0.0012 
0.0014 


0.0014 
0.0016 








Term Estimate | Std Error t|t Ratio Prob>|t| 
Ratio 


(TwrClassRng-123.651)*(CommCentRel- -0.000112] 3.489e-5 — 0.0019 
66.4419) 


(MiniUAVFuel-2016.79)*(TwrHt-34.8217) |0.0000159| 4.975e-6| 3.19/- 1m |] 0.0020 
(MiniUAVFuel-2016.79)*(PredPClass-0.55372) _[ -0.000647| 0.000207| -3.13/- maj 1] | 0.0025 
(CommCentLat-10.2403)*(NumBPCars-3.51938) | -0.013412] 0.004459[ -3.01|[ eT | (0.0035 
TwrHt-34.8217)*(TwrPClass-0.55698 | -0.035852| 0.0125] -2.87/[7 dC. 


(MiniUAVFuel-2016.79)*(MiniUAVCommLat- ee | 0.0062 
MInlUAV Turret 76.651) CommCentFer [po ites | kes 26 ——— 0.0077 
IMU TurvetF-O75 651) WAnTUAVCommPer pees eee 8 ee 0.0079 
ProaciaasFing- 134.500)" (CommCanter $<" 0.0122 
IMU TurvetF-S75 651)" MInUAVGlassAng. | — 0.0259 
IMintUAVGlassBng 1.8682)*(PredClassRno- 0.000268 0.000115) -2.24 0.0277 
emcees 23.651)*(CommCentLat-  P.00242/ 0.000111) -2.18 [| 0.0320 
Mina SALTO TERY (CommCentRel- ee SS  K§¥——===—| 0.0332 


(MiniUAVFuel-2016.79)*(TwrPClass-0.55698) -0.000474| 0.000226] -2.10/- ms __]]| 0.0391 








Table 16. Linear Regression Results for Total Number of 
Classified Immigrants with MiniUAV 





4. Total Number of Classified Smugglers 
a. Without MiniUAV 


Although the p-value is less than 0.0001, 





regressors are less successful at explaining the overall 
variance of the MOE, with an RSquare value of 0.74, than 


the model for the number of classified immigrants. 








Despite the lower RSquare value, the important 


factors are in accordance with the model for the number of 





classified immigrants. Tower classification range becomes 








the most important term, instead of the Predator 


classification range, and the number of BP cars emerges as 





an important factor for this MOE (Table 17). 
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Response Mean(TotSmugClassified) 
Actual by Predicted Plot 





Mean(TotSmugClassified) 
Actual 





1.6 1.7 1.8 1.9 2.0 


Mean(TotSmugClassified) Predicted 
P<.0001 RSq=0.74 RMSE=0.0433 








Summary of Fit 


RSquare 0.7413 

RSquare Adj 0.716768 

Root Mean Square Error 0.043339 

Mean of Response 1.933854 

Observations (or Sum Wgts) 128 

Analysis of Variance 

Source DF Sum of Squares Mean Square F Ratio 
Model 11 0.62431321 0.056756 30.2177 
Error 116 0.21787429 0.001878 Prob > F 
C. Total 127 0.84218750 <.0001 


Sorted Parameter Estimates 





Term Estimate Std t|t Ratio Prob>|t| 
Error| Ratio 

TwrClassRng [0.0010978| 8.863e-5| 12.39|[ | <.0001 

PredClassRng 10.0008979| 8.939e-5| 10.05/[- | I] <.0001 

(TwrClassRng-126.367)*(PredClassRng-125.805) | -1.483e-5| 2.223e-6| -6.67|[ mm 7 «dC <.0001 

NumBPCars [0.0128331| 0.004003] 3.21|[ Tm id SC(0.0017 


(PredClassRng-125.805)*(CommCentRel- -1.18e-5] 4.802e-6 | 246 0.0154 
63.5547) 


(CommCentRel-63.5547)*(NumBPCars-3.48438) | -0.000504| 0.000209] -2.41|— —____jm Sf ____]| 0.0175 








TwrPClass 0.043929| 0.018993] 2.31/[[ TT rm —idC(0.0225 

TwrHt 0.0009921| 0.000443] 2.24/T ST | _0.0270 

(TwrHt-35.1875)*(PredClassRng- 125.805) -2.541e-5| 1.141e-5] -223/- mm 1 _]| 0.0279 
Table 17. Linear Regression Results for Total Number of 


Classified Smugglers without MiniUAV 


b. With MiniUAV 


The model for the number of classified smugglers 
with MiniUAV explains less variability in the MOE, with an 


Rsuared value of 0.88, compared to the model for the number 





of classified immigrants. However, important factors are in 
719 








accordance with the without miniUAV scenarios. The main 








difference between these two cases is the order of the 








first three most important factors. Tower classification 


range becomes the most important term instead of Predator 





classification range (Table 18). 





Response Mean(TotSmugClassified) 
Actual by Predicted Plot 





1.9 


1.7 


Mean(TotSmugClassified) 
Actual 
b 


=e 
oO 





1.60 1.70 1.80 1.90 2.00 


Mean(TotSmugClassified) Predicted 
P<.0001 RSq=0.88 RMSE=0.0367 








Summary of Fit 


RSquare 0.88387 
RSquare Adj 0.831083 
Root Mean Square Error 0.036707 
Mean of Response 1.934109 
Observations (or Sum Wgts) 129 


Analysis of Variance 


Source DF Sum of Squares Mean Square F Ratio 
Model 40 0.9024606 0.022562 16.7442 
Error 88 0.1185730 0.001347 Prob > F 
C. Total 128 1.0210336 <.0001 


Sorted Parameter Estimates 





Term Estimate Std t|t Ratio Prob>|t| 
Error| Ratio 

TwrClassRng 10.0013447| 0.000116/ 11.64/- _____f_ | <.0001 
(TwrClassRng-123.651)*(PredClassRng-124.209) | -2.246e-5| 2.593e-6| -8.66|[ [ummm 1 || <.0001 
PredClassRng |0.0007661| 0.000101] 7.62|[ | mmm | <.0001 
(MiniUAVCommLat-10.1783)*(CommCentRel- 0.000212) S-7S8e-5) “6.67 | Ty <.0001 
66.4419) 

(MiniUAVFuel-2016.79)*(TwrHt-34.821 7) | 3.81e-6| 8.577e-7| 4.44/[—" | mm ____]| <.0001 


(MiniUAVCommRel-79.3256)*(CommCentLat- -0.000335| 7.649e-5 —— <.0001 
10.2403) 
(MiniUAVTurretHt-275.651)*(NumBPCars- 0.0002075| 5.23e-5 a 0.0001 
3.51938) 


(MiniUAVFuel-2016.79)*(TwrClassRng-123.651) | -5.161e-7 LS2te7| -3.39{T 0.0010 


(MiniUAVCommLat-10.1783)*(TwrClassRng- 6.0182e-5| 1.78e-5 ee 0.0011 
123.651) 


(MiniUAVTurretHt-275.651)*(TwrHt-34.821 7) -2.023e-5| 6.134e-6] -3.30|/" mmm 1 | 0.0014 
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Term Estimate Std t|t Ratio Prob>|t| 
Error| Ratio 


(MiniUAVTurretHt-275.651)*(CommCentRel- -0.000011] 3.372e-6 2 0.0014 
66.4419) 


(TwrHt-34.8217)*(CommCentRel-66.4419) 0.0000926| 2.94e-5[ 3.15|/- 1 mmm _]| 0.0022 


(MiniUAVClassRng-11.8682)*(PredClassRng- -0.00006] 1.946e-5 SS 0.0028 
124.209) 
(MiniUAVTurretHt-275.651)*(TwrClassRng- -4.708e-6| 1.553e-6 a 0.0032 
123.651) 


TwrPClass 0.0682418| 0.0236] 2.89/- 1 mmm __]| 0.0048 
MiniUAVTurretHt 0.0001615] 5.873e-5|  2.75[—> 1 mm | 0.0072 
(MiniUAVClassRng-11.8682)*(CommCentLat- 0.0004365| 0.00016 — 0.0078 
10.2403) 

MiniUAVFuel 1.7017e-5| 6.274e-6| 2.71/- 1m | _0.0080 


(MiniUAVFuel-2016.79)*(PredPClass-0.55372) -6.576e-5| 0.000025| -2.62/[ Sm 7 —C~dY:«C«O08 
(MiniUAVFuel-2016.79)*(PredClassRng-124.209) | -3.064e-7[ 1.172e-7| -2.61/~ —___mm 1 | 0.0105 


(MiniUAVFuel-2016.79)*(MiniUAVCommRel- | 0.0158 
(Sa CeniLaL T0205" Commenter Pee | 0.0194 
IMnIUAV Turret 75 651 )*(CommCentLat- | -$.075e5] 181665) 2.94 0.0218 
TMiniUAVGommAel-75 S256)" ProdPClaes: 0.004284) 0.001881) 2.91 0.0231 
(AIA SRT SEY TCR | ——————— 0.0231 








CommCentLat -0.001787[ 0.000799 -2.24[- my 7 _]| 0.0278 

(PredPClass-0.55372)*(CommCentRel-66.4419) | -0.001692] 0.000764] -2.21/- mm 1 ___]| 0.0293 

CommCeniRel 0.0003503|0.000171[_ 2.05[- 1 mm _]| 0.0435 
Table 18. Linear Regression Results for Total Number of 


Classified Smugglers with MiniUAV 


D. REGRESSION TREES 


Regression, or partition, trees are another way of 





analyzing classification and regression problems. Because 
regression trees are constructed by if-then conditions, it 
is non-parametric and interactions between the factors 
emerge automatically. Regression trees present some 


threshold values which are hard to determine by linear 








regression, and also it is easy to interpret the results in 


most cases. 





Regression trees for all MOES are presented in 


Appendix B. 
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1. Total Number of Captured Illegal Entrants 
a. Without MiniUAV 


The important factors emerged from linear 
regression, such as the number of BP cars, communication 


center reliability and tower classification range, were in 





accordance for the number of captured immigrants and 
smugglers. The result from partition trees for these 
scenarios complies with linear regression. Because number 


of BP cars is too dominant for these MOEs, the tree 





initially splits on all possible levels of this factor. The 
mean of the captured illegal entrants for the scenarios 


with more than three BP cars appears to be larger than the 





mean for the scenarios with fewer BP cars. While the 
threshold for communication center reliability for 
immigrants is around 40 for the scenarios with three and 


four BP cars, it is 73 for the scenarios with five BP cars. 














This indicates that more reliable communication systems are 





needed when the manpower is increased. 


b. With MiniUAV 


Although the number of BP cars is the most 
important factor for the scenarios with the MiniUAV too, 


some additional factors also are very important and come 





into play for these scenarios, such as MiniUAV 





classification range and MiniUAV endurance. Also, it is 





possible to see some differences between the trees for 


immigrants and smugglers. While the endurance of the 





miniUAV seems to be as important as miniUAV classification 





range for capturing immigrants, classification parameters 





for the sensors (miniUAV, towers and Predator B) become 


more important for capturing smugglers. The threshold value 
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for miniUAV endurance is roughly 90 minute (minimum of 60 





minutes, maximum of 120 minutes). It becomes relevant to 





hav xtended endurance for the miniUAV when its 








classification range is more than 500 meters, which means 


that endurance and classification range should be improved 





simultaneously. 
2. Total Number of Classified Illegal Entrants 
a. Without MiniUAV 


The most important factors emerging from the 
partition tree for the total number of classified 


immigrants and smugglers are Predator B classification 











range and tower classification range, just as the results 


from linear regression. With their large classification 








ranges, it is normal to see these factors in the partition 


tree; however, they have unexpected interaction effects. If 











the classification range for the Predator B is low, the 


classification range for tower needs to be high. However, 





if the classification range for the Predator B is high, the 











classification range for the tower becomes less important. 





The same relationship is true for Predator B, as with the 





cases of tower classification range. This result is an 





indication of negative correlation between the 


classification ranges of the Predator B and towers. This 





negative correlation can be seen by looking at the negative 





coefficient of the interaction term of these factors in the 


linear regression models. Due to the strong relationship 





between these two factors, it is not possible to determine 


stand-alone threshold values. However, 5000-6000 m for the 





towers and 5000-6500 m for the Predator B can be accepted 
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as the threshold values that make significant differences 


for the mean of classified illegal entrants. 





b. With MiniUAV 


Similar to the scenarios without the miniUAV, the 








classification ranges of towers and Predator B are the most 
important factors and they have a strong correlation. The 


relatively small classification range of the miniUAV, 





compared to the Predator B and towers, can explain the 
reason that the miniUAV does not emerge as an important 


factor. The threshold values for the classification range 








of the towers are 3300 m for classified immigrants and 5500 





meters for classified smugglers. 
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VI. CONCLUSION 


To achieve an improved border security, the CBP 
focuses on improving manpower, technology and 


infrastructure along the border. Due to the long land 





border and limited assets, it is difficult or impossible to 





fully protect every point along the border. Therefore, CBP 
must direct its resources in an intelligent way, 
determining the assets to invest in and the best way to 


employ them. 


A border model is constructed in MANA for this study 





to explore the effectiveness of the assets used in border 
security. According to the analysis, three important 
factors or systems need to be taken into account for an 
effective border security system: Manpower, MiniUAV usage 


and communication systems. 


An infrastructure such as a fence along the border or 


a state of the art sensor system may be used to prevent 





illegal entrants. However, although it is very expensive to 
construct an infrastructure, it is easy to destroy it or 


find a way to avoid it if there is nothing to protect this 








infrastructure. Similarly, illegal entrants may not care if 





they are detected by the sensors when there is no physical 








deterrence to stop them. The only asset that can protect 


the infrastructure or provide a physical deterrence is 





manpower, particularly BP agents. Without a sufficient 
number of BP agents to assist the overall security system, 


other measures becom ineffectiv in stopping illegal 








activities. Particularly for this model where 50 km of land 
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border is used, an average of seven BP agents can only 





capture roughly 25% of the 22 illegal entrants. 


Comparison tests for various MOES reveal that using a 


hand-launched miniUAV improves th ffectiveness of the 





system. Because miniUAVs are controlled directly by BP 





agents, they need less coordination effort. Any other 
sensor that can be used directly by BP agents may give 
Similar results, but airborne sensors have some advantages 
compared with ground sensors. For example, it is not 


necessary to acquire land to construct any facility, and 





this fact is especially important in the modeled region, 
which contains many sensitive lands. Also they can track 


the illegal immigrants without any need for a road in rough 








terrains. To improve the border security system, MiniUAV 








classification range and endurance ar the two most 





important factors that emerge in the regression models. For 


the sensors other than the miniUAV, classification ranges 








of the Towers and Predator B are other important factors 


regarding the border security system. 


Having found the usage of the miniUAV to be an 
important contribution to the border security system, 


further analysis should be done by developing more specific 





models to reveal th effectiveness of different tactics, 


techniques, and procedures (TTPs). 


A reliable communication web is another important 





factor regarding an effective border security system. 





Having perfect sensors that can detect and classify illegal 
activities is not sufficient if this information is not 
successfully transmitted to the BP agents. Better 


communication systems should be sought, which integrate 
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sensors with each other and the BP. Other than reliability, 





communication latency, which is the time to transmit data 





from one agent to another, is another factor that emerges 
as important in the linear regression models. Further 


analysis can be done to find the critical levels of these 





important factors by using higher-resolution models and 





examining their marginal contribution to the MOEs. 


Border protection is a vital national security issue. 





Although gaps in the security system of the border 








contribute to the surge in the number of illegal activities 


along the border, it is not possible to completely stop all 








illegal trafficking, even with a perfect security system. 








But this fact does not mean that the number of the illegal 








activities will remain the same regardless of the security 





measures that are taken. A better security system decreases 
the number of the illegal activities by deterring the 
illegal entrants. And this thesis has explored the use of 








UAVs to enhance border security with the results indicating 





the potential effectiveness of these assets. 
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APPENDIX A: EXPERIMENTAL DESIGN 


proaasac) | | | | 
Z| <|2Z|<|Z 
8}°|81° |withUAV?(Y/N) 


o 


an 
= 
i aa 


Zz 
°o 


= 
NI 
De) 
ice) 
o 


alves None [No 
i[Yes [Yes [No_[Yes | 
— 4 0.32 ; ; | 19] 4 OlYes [No_[No_[No_| 

1 93 -97{_ 81] 54] 8 iYes [Yes [No_[Yes | 
Tro S02 1a 19 97 83] $ 92 .96] 47] 31 ofYes [No_[No_[No_| 
5 i[Yes [Yes [Yes [Yes | 

Q[No_[No_[No_[No_| 
iYes [Yes [No_[Yes | 
49 
[ 0.97] 7 6 . i[Yes [Yes [No_[Yes | 

rats oe[ a ee af ze) oat asl ef sal rsa oars oa ros ral oe 28-17 O[No_[No_[No_[No_| 


=|/o]~NI 
QO] CO} 
=|) 
E 
« N 
= 
=I 
rS 
iY) 
=| a] 
nN) 


=] 
DID 







































































13 1 
east iat 7 taal atonal os ol raad oa a[asa 0 
i 
0 

Paral [0.76 . 3] 0 
a al-rataa[-oa[ fea] 207|- 0-8] Teal Sot 0.68 . . 5 0 
Fa Ee [— 29 1 

97 


O[Yes [No [No | 
0 
1 


rasali9 roa EE 
Bom 
Bal aat| Tel 2a] 20 
EOEIEL toasts : 58 T[Yes |No_|[No_|Yes | 
[338] 8 12] 19] [| 0.51] 17: ; En 94] 36[ 24] 16] 71 OlYes [Yes [No_[No_| 
sor[ore] 91a 1 [ 29 ai E Ba i] 0. ; 65|_1| 52 T[Yes [No_[No_|Yes | 
: | oles [Yes [No_[No_| 
1[Yes [No_[No_[Yes | 
olYes [Yes [No_[No_] 
. 0.26 . ; T[No_[No_[No_[Yes | 
| 0.3] 82[ 123] 0.59 7{ 41[ 27 3[ 53 i[Yes [No_[No_[Yes | 
C 





o@]—lo 
@/@]—]|r 
ALOLO|H 

als 

alt 


No 


&IS| SI 
ce) Kes) Ie) 
| 
hyo} — 
an 


= 
ie.) 
al 


Bla 00] 
SEE 


1 


oO 

BERBER UUERRUBEE 

CIDIO[DIOJDIOJOJO/MIOJoOJO/VIoO]aIo 
oO Oo Oo oO oO oO 


x 


|= 
ole 
ARLOLO 
a] 



































o 
£ 
5 
S 
g | | 
we 
— 
-— 
-— 
-— 
-— 
is 
—4 
-— 
-— 
-— 
-— 
-— 
-— 
Ll 
-— 
-— 
-— 
-— 
-— 
-— 
7 
-— 
-— 
-— 
-— 
-— 
-— 
-— 
-— 
=a 








95 





Factor Name 
IMiniUAV Fuel 
IMiniUAVDetRng 
iTwrClassRng 
iTwrPClass(2) 
TwrDetRng(2) 
iTwrClassRng(2) 
iTwrPClass(1) 
TwrDetRng(1) 
iTwrClassRng(1) 
redClassRng 
PredDetRng 
PredPClass(2) 
PredDetRng(2) 
PredClassRng(2) 
PredPClass(1) 
PredDetRng(1) 
PredClassRng(1) 
IWithUAV? 
BPCar4(Sqd8) 


Z| 
oO 
Z| 
°° 
Z| 
° 


94] 37] 0.82 . - K 19 
86] 49] 0.86] 6: . - a 58 
82] 45] 0.68 





colon 
= 





a 
8 
6 


1213 
123 
1381 


to} 
ES 
s|°o 








171 
281 
221 








alo/a 


Sake) 


8] 0.59 
3] 36] 0.8 
8] 34] 0.66 














0.94 


171 | E 4 R .55 0.89 
1494] 246} 10 16 | 4 i : 

36] 0.66 78| 0.86] 5 

32] 0.32 189} 0.62] 126 








oO} 
als 
glo 
als 
gle 
<|Z 
D|Oo 
a 





~<| 
fo) 
Oo 
=< 
fo) 
o 
Z| 
° 
a 








s]o|/-|0 








0.4 
0.47 














1 0.92 


8 
95} 63] 0.95 


| 80] 0.84] 60] 40] 2 


59] 0.97 
114] 0.91 
73] 0.85 
155] 0. E Yes |No_|No_|No 
209] 0. E [Yes [Yes [Yes [Yes 
159] 0. H Yes |No_|No_|No 








~<| 
o 
Oo 
~<| 
fo} 
Oo 
Z| 
°° 
A 


Yes |No_|No 
Yes |Yes |Yes 
No _|No_|No 








o|-[0 








1794] 230] 19] 29] 17 

913} 322] 11] 17| 18 

1550} 162] 9] 14] 19} 75 
573 





=o] =| 





+|° 
~<| 
fo) 
oO 
~<| 
fo} 
oO 
Z| 
° 
a 


o|= 








Yes |No 
Yes 





0.96 
0.6 





o|+|0o 
Z| 

o 

Zz 

o 


Yes 
Yes |Yes 


fo) 


0.89 
0.92 
0.8 


i=} 
Z| 
°° 
Z| 
°o 





Yes 
Yes |No 


Yes |No_|No_|No 

Yes |Yes [Yes |Yes 

Yes |No_|No_|No 
[Yes | 


=< 
) 
Oo 





s]° 








o|-|0 


=< 
) 


Yes |Yes 


135] 9! 0.9 No_|No 
0.85} 89] 59] 0.92 Yes |Yes 








0.66] 182] 121] 0.76 




































































1513] 199 
1053] 357| 16 
1213] 176] 17 


39 
138 


Factor Name 
IMiniUAV Fuel 
MiniUAVTurretHt 
IMiniUAVClassRng 
IMiniUAVDetRng 
MiniUAVCommLat 
ITwrClassRng 
rDetRng 
iTwrPClass(2) 
iTwrDetRng(2) 
iTwrClassRng(2) 
iTwrPClass(1) 
TwrDetRng(1) 
iTwrClassRng(1) 
redClassRng 
PredDetRng 
PredPClass(2) 
PredDetRng(2) 
PredClassRng(2) 
redPClass(1) 
redDetRng(1) 
redClassRng(1) 
IWithUAV? 
BPCar4(Sqd8) 
WithUAV?(Y/N) 


Qao;a 








NI 


y for) eal 
to] 
|] 00 

= 
tol NB] ola 


an 

G 

(Ss) 

a 

a 

f= 

Ss 

Zz 
5 Yes 
2 No 
4 Yes 


149 Yes 
107] 0. : i Yes 
206] 0. ; | Yes 
171} 0. 
281 E 1 Yes 
221 Yes 
Yes 
No 
Yes 
i|Yes 
Yes N 
No N 
Yes N 
Yes |No_|N 
Yi 

lo_|N 

Ni 























[e) 
[e) 
10 
[e) 





Yes s [No 


e 
0|Yes lo 
il 





No 

Yes 
Yes 
Yes 
Yes 


O|Dlo 
Oo 





Yes 
i|Yes 





182 I Yes 
249 Yes 
201 i | Yes 
2 Yes 
0.9 Yes 
A Yes 
Yes 














d 


) 
oa 








=<]<[~< 


i Yes 
0.69] 104] 69] 0.79 Yes No 
i|Yes Yes 
Yes No 
Yes No 
20} 80} 0.83 | No No 
Yes No 


























Yes 
Yes No 
i|Yes No 
No 
Yes 











|S 
oO 


Yes 
No 

Yes 
Yes 
Yes 
Yes 
Yes 





a 
o 


=| 
oO 
8 





° 





</Z 
fo) 
a 








0.7[173[115] 0.8 
0.68] 54[ 36] 0.78 


101 
99 
83 


4 R 56 i 89 
46] 0.29 7 77 i 0.66] 182] 121 
35] 0.35} ql 92 B 0.72] 89) 5 
















































































redPClass 
redClassRng 


Factor Name 
IMiniUAV Fuel 
MiniUAVTurretHt 
IMiniUAVClassRng 
IMiniUAVDetRng 
MiniUAVCommRel 
jTwrHt 

‘wrPClass 
ITwrClassRng 
iTwrDetRng 
iTwrPClass(2) 
iTwrDetRng(2) 

wrClassRng(2) 
iTwrPClass(1) 
TwrDetRng(1) 
iTwrClassRng(1) 
PredDetRng 
PredPClass(2) 
PredDetRng(2) 
PredClassRng(2) 
PredPClass(1) 
PredDetRng(1) 
PredClassRng(1) 
IWithUAV? 
BPCar5(Sqd9) 
BPCar1(Sqd7) 
BPCar4(Sqd8) 


IT 
Pi 
Pi 











0.66 
0.8 











20 
21| 0.57 
2| 0.77 








0} 0.91 
174] 116] 0.86 
uA 
1 i 























s]o|/+|o 








o|-[0 








251| 0.92] 167] 111 
252| 0.91] 168] 112] 0.96 
128] 0.74] 86] 5 

288 0.7 
50} 0.93 
48] 0.95 
9] 0.77 














2619 18] 27 
2834] 152) 11] 17 
3181] 362] 10 


0.9 
0.97 


34] 0.3 k Yes 
72| 36} 0.44 i Yes 
77| 48] 0.43 } E | Yes 


Yes 
Yes 
Yes 




















rs 
oO 
ra 
oO 





0.83 } Yes 

3 Yes 
Yes 
Yes 
Yes 





<|Z| 
2|o 








o|+|o 
Z| 
°o 
Z2|Z 
fou Ke) 
Z2\Z2 
o|o 


I Yes 
195} 0.83 0.91 5 - Yes 
258] 0.83 0.91 1 0.51} 129) 8 

243] 0.93 0.97 0.91] 84] 56] 0.96 
254] 0.9 0.96 0.91} 116] 77] 0.96 


2928] 15 | i 62] 41 
3191] 372] 13 75) 191] 127] 0.95 
2 | 5 18 | 5 128] 85] 0.95 


277: 


2188 f 
2497} 295] 20 14] 82] 47] 0.25 
3 74 i 


fo) 


5 
rs 
oO 
ia 
oO 








<| 
fo} 
ao 
Z2)/Z)2) 
O|O|o 








=< 
() 
a 








A 
°o 






























































2|Z 
o|o 








IMiniUAVDetRng 
MiniUAVCommLat 
‘wrPClass(2) 
wrDetRng(2) 
wrClassRng(2) 
iTwrPClass(1) 
PredClassRng 


c 

=| 9 
S| 
x| 3 
o/s 
Q/o 
=| 
E 


Factor Name 
PredPClass 
PredDetRng 


|MiniUAVCommRel 
iTwrDetRng 


3 |MiniUAVFuel 
Q|TwrHt 





= 
nN 
o 
oO 





PredPClass(2) 
PredDetRng(2) 
PredClassRng(2) 





0.5 


PredPClass(1) 
PredDetRng(1) 
PredClassRng(1) 
(CommCentLat 
CommCentRel 
INumBPCars 

o|WithUAV? 


for) 


BPCar5(Sqd9) 





| B|MiniUAVTurretHt 


wn 








0.88 











[155] 1 
3 





131 
131 


[103] 
87 
87 


1|Yes 
OlYes [No_|No 


-4[__1\¥es [Yes |No_|ves 


|_O|No_[No_[No_| 





jo 
1|Yes |Yes |Yes 

O]Yes [No_|No 
S 






























































\WithUAV?(Y/N) 






<]2Z|<|/Z 


ala (4 


THIS PAGE INTENTIONALLY LEFT BLANK 


100 





APPENDIX B: REGRESSION TREES 








































































































































































































A. Regression Tree for the Mean Number of Captured 
Immigrants 
RSquare N Number of 
Splits 
0.916 128 8 
All Rows 
Count 128 LogWorth Difference 
Mean 3.0304687 82.535553 1.84046 
Std Dev 1.0914607 
NumBPCars<4 NumBPCars>=4 
Count 65 LogWorth Difference Count 63 LogWorth Difference 
Mean 2.1246154 31.888939 0.93961 Mean 3.9650794 10.584178 0.84524 
StdDev 0.5417076 Std Dev 0.6240528 
NumBPCars<3 NumBPCars>=3 NumBPCars<5 NumBPCars>=5 
Count 22 || Count 43 LogWorth Difference Count 42 LogWorth Difference Count 21 LogWorth Difference 
Mean 1.5030303 |} Mean 2.4426357 —2.5485098 0.44339, Mean 3.6833333  3.4558445 0.57167 Mean 4.5285714 —1.1996543 0.55962 
Std Dev 0.2462221 || Std Dev 0.3330416 StdDev 0.478664 Std Dev 0.4873153 
CommCentRel<41 ||CommCentRel>=41 CommCentRel<42 CommCentRel>=42 CommCentRel<73_|/CommCentRel>=73 
Count 7 || Count 36 Count 12 Count 30 LogWorth Difference Count 13) Count 8 
Mean 2.0714286 || Mean 2.5148148 Mean 3.275 Mean 3.8466667 1.0004487 0.4236 Mean 4.3153846 |) Mean 4.875 
Std Dev 0.2189441 || Std Dev 0.3032569 Std Dev 0.3908357 Std Dev 0.4115758 Std Dev 0.3207674|| Std Dev 0.529675 
TwrPClass>=0.35 TwrPClass<0.35 
Count 


23 LogWorth Difference || Count 
Mean 3.7478261 1.3064686 


Std Dev 0.3692507 


7 
0.43137 |] Mean 4.1714286 


StdDev 0.397878 





























TwrPClass<0.68 TwrPClass>=0.68 


Count 17 || Count 
Mean 3.6352941 || Mean 
Std Dev 0.3344711 


6 
4.0666667 
Std Dev 0.2772885 
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B. Regression Tree for the Mean Number of Captured 


Immigrants with MiniUAV 






















































































RSquare N Number of 
Splits 
All Rows 
Count 129 LogWorth Difference 
Mean 5.4937984 63.779265 3.06241 
Std Dev 1.89234 
Num BPCars<4 Num BPCars>=4 
Count 63 LogWorth Difference Count 66 LogWorth Difference 
Mean 3.9269841 23.686162 1.61905 Mean 6.9893939 14.760045 1.73373 
Std Dev 0.9754304 Std Dev 1.2221203 
Num BPCars<3 Num BPCars>=3 MiniUAVClassRng<10 MiniUAVClassRng>=10 
Count 21 LogWorth Difference Count 42 LogWorth Difference Count 24 LogWorth Difference || Count 42 LogWorth Difference 
Mean 2.847619 2.9064608 0.81 Mean 44666667 8.8487607 0.89286 Mean 5.8861111 2.8661519 0.76476 || Mean 7.6198413 5.4892153 1.3 
Std Dev 0.5842917 Std Dev 0.6141511 Std Dev 0.6110035 Std Dev 1.0191941 



































MiniUAVFuel<1813 |/MiniUAVFuel>=1813 
Count 10 |] Count 1 
Mean 2.4233333 || Mean 3.2333333 











MiniUAVFuel<1747 MiniUAVFuel>=1747 MiniUAVFuel<1663 |/MiniUAVFuel>=1663 
Count 14] | Count 28 LogWorth Difference Count 10 








MiniUAVFuel<2384 ||MiniUAVFuel>=2384 









































Count 14 Count 27 || Count 15 
Mean 3.8714286 | | Mean 4.7642857 4.246058 0.6975 Mean 5.44 || Mean 6.2047619 Mean 7.1555556 || Mean 8.4555556 
Std Dev 0.4594817 || Std Dev 0.3861491 Std Dev 0.3393954| | StdDev 0.4916528 Std Dev 0.4577994 || Std Dev 0.5029583 Std Dev 0.8192659]| Std Dev 0.7946914 


























MiniUAVClassRng<16]| MiniUAVClassRng>=16 
Count 20 Count 8 
Mean 4.565 Mean 5.2625 
Std Dev 0.3661839 Std Dev 0.4138236 
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C. 


Smugglers 








All Rows 

Count 128 
Mean 0,2984375 
Std Dev 0,1021946 





LogWorth Difference 
7,7160486 0,08956 











RSquare N Number of 
Splits 
0,489 128 9 
NumBPCars<4 
Count 





Std Dev 0,0904405 


65 LogWorth Difference 
Mean 0,254359 3,4747704 


0,07759 























Num BPCars<3 
Count 
Mean 


Std Dev 0,0748254 


0, 


22 LogWorth Difference 


2030303 3,7490092 








Num BPCars>=3 

Count 43 
0,10142 Mean 0,2806202 
Std Dev 0,0870542 























CommCentLat<10 


Count 


Mean 0,1615385 
Std Dev 0,0487741 


13. LogWorth Difference 
1,0935731 0,0627 








CommCentLat>=10 
Count 9 
Mean 0,262963 
StdDev 0,065499 




















Count 





CommCentRel<57 


Mean 0,1277778 
Std Dev 0,0443053 


6 


CommCentRel>=57 
Count 7 


Mean 0,1904762 








Std Dev 0,0317063 











Num BPCars>=4 
Count 

Mean 0,3439153 
Std Dev 0,0938666 





63 LogWorth Difference 


1,6504561 0,05556 











Regression Tree for the Mean Number of Captured 

































































NumBPCars>=5 

Count val 
Mean 0,3809524 
Std Dev 0,0840446 








NumBPCars<5 

Count 42 LogWorth Difference 

Mean 0,3253968 | 1,3585154 0,08611 

StdDev 0,0939381 
TwrPClass<0,39 TwrPClass>=0,39 
Count 12 LogWorth Difference Count 30 LogWorth Difference 
Mean 0,2638889 1,841644 0,10762 Mean 0,35 2,5123205 0,10181 
Std Dev 0,0717154 Std Dev 0,0912871 
PredClassRng>=110]|PredClassRng<110 TwrHt<37 TwrHt>=37 
Count 7 |] Count 5] | Count 17 LogWorth Difference || Count 13 
Mean 0,2190476 || Mean 0,3266667 


Std Dev 0,0465759 











StdDev 0,0494413 








Mean 0,3058824 
StdDev 0,0637217 


0,8088503 0,07095 








Mean 0,4076923 
StdDev 0,091443 
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CommCentRel>=48 
Count 10 
Mean 0,2766667 
Std Dev 0,0609594 





CommCentRel<48 
Count 7 


Mean 0,347619 








Std Dev 0,0424139 














D. 


RSquare 


0.589 


Regression Tree for 


Smugglers MiniUAV 


129 


Number of 
Splits 
9 


the 





All Rows 
Count 
Mean 





129 LogWorth Difference 
0.5129199 6.0894644 
Std Dev 0.1566971 


0.16526 

















Count 
Mean 





Num BPCars<5 


Std Dev 0.1388241 


107 LogWorth Difference 
0.4847352 4.4127929 0.11177 














Count 





MiniUAVClassRng<11 


46 LogWorth Difference 














MiniUAVClassRng>=11 











Count 61 LogWorth Difference 
Mean 0.4210145 = 1.9120072 0.08962 Mean 0.5327869 6.9110967 0.19897 
Std Dev 0.1209752 StdDev 0.1327057 
Num BPCars<4 NumBPCars>=4 MiniUAVFuel<1250 
Count 














Mean Number of Captured 





NumBPCars>=5 
Count 

Mean 0.65 
Std Dev 0.1692655 


22 LogWorth Difference 
1.7114813 0.18485 























MiniUAV Class Rng<12||MiniUAVClassRng>=12 
Count 1 Count 11 LogWorth Difference 
Mean 0.5575758 Mean 


Std Dev 0.1022969 








0.7424242 
Std Dev 0.1758098 


0.8416511 0.22778 











26 LogWorth Difference 
Mean 0.3820513 — 3.0400417 














Count 20 Count 














MiniUAVFuel>=1250 


























11 || Count 50 LogWorth Difference 
0.14206 Mean _0.4716667| | Mean 0.369697 || Mean 0.5686667 _2.0379421 0.11151 
Std Dev 0.1159428 Std Dev 0.1104351 Std Dev 0.0912318 || Std Dev 0.1121648 
TwrClassRng<121 || TwrClassRng>=121 Num BPCars<3 Num BPCars>=3 
Count 12 || Count 14 LogWorth Difference Count 8 || Count 42 
Mean 0.3055556 |} Mean 0.447619 1.3744468 0.10519 Mean 0.475 |] Mean 0.5865079 
Std Dev 0.0885727 || Std Dev 0.0958224 Std Dev 0.1433721 || StdDev 0.097408 
Num BPCars<3 Num BPCars>=3 
Count 5 | Count 9 
Mean 0.38 |} Mean 0.4851852 
Std Dev 0.0730297 || Std Dev 0.0883665 
































104 









































PredPClass<0.62 
Count 





Mean 0.6388889 |} Mean 
StdDev 0.118165 |] Std Dev 


PredPClass>=0.62 
6 || Count 5 
0.8666667 
0.1563472 

















E. Regression Tree for the Mean Number of Immigrants 


Classified 











ee ea eee 





















































RSquare N Number of 
Splits 
0,863 128 7 
All Rows 
Count 128 LogWorth Difference 
Mean 18,905469 19,364874 1,20349 
Std Dev 0,9787039 
PredClassRng<136 PredClassRng>=136 
Count 74 LogWorth Difference Count 54 LogWorth Difference 
Mean 18,397748 31,515726 1,56072 Mean 19,601235 8,0362699 0,35379 
Std Dev 0,9983294 Std Dev 0,2633411 
TwrClassRng<127 TwrClassRng>=127 TwrClassRng<105 |/TwrClassRng>=105 
Count 36 LogWorth Difference Count 38 LogWorth Difference Count 17 || Count 37 LogWorth Difference 
Mean 17,596296 6,3108174 1,27024 Mean 19,157018 6,023368 0,49185 Mean 19,358824 || Mean 19,712613  1,6923712 0,14353 
Std Dev 0,7989985 Std Dev 0,3813921 Std Dev 0,2916737 || Std Dev 0,1542307 


























PredClassRng<73_ ||PredClassRng>=73 PredClassRng<92_ |/PredClassRng>=92 PredClassRng<174 || PredClassRng>=174 

Count 8 || Count 28 Count 18 || Count 20 LogWorth Difference Count 20 || Count 17 
Mean 18,898148 || Mean 19,39 4,1432628 0,44889 Mean 19,646667 || Mean 19,790196 

Std Dev 0,3278913 || Std Dev 0,2577682 Std Dev 0,1690704 |} Std Dev 0,0880062 


Mean 16,608333 |} Mean 17,878571 
Std Dev 0,6311365 || Std Dev 0,5936713 




















TwrClassRng<147 ||/TwrClassRng>=147 
Count 5 || Count 15 
Mean 19,053333 |} Mean 19,502222 
Std Dev 0,2340465 || Std Dev 0,1433536 
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F. 


Regression Tree for the Mean 


Classified with MiniUAV 


RSquare N 


0,821 129 


Number of 
Splits 
6 





All Rows 
Count 
Mean 


129 


19,020672 17,453178 


Std Dev 0,9481634 


LogWorth Difference 


1,73836 











TwrClassRng<66 




























Number of 























Immigrants 












































TwrClassRng>=66 
Count 15 LogWorth Difference Count 114 LogWorth Difference 
Mean 17,484444 4,2617229 2,24259 Mean 19,222807 18,459012 0,82744 
Std Dev 1,3592053 Std Dev 0,6595239 
a 
PredClassRng<100 ||PredClassRng>=100| |PredClassRng<119 PredClassRng>=119 
Count 6 || Count 9 Count 49 LogWorth Difference Count 65 LogWorth Difference 
Mean 16,138889 |} Mean 18,381481 Mean 18,75102 12,486933 1,00275 Mean 19,578462  15,272957 0,45467 
Std Dev 0,9684849 || Std Dev 0,6193884 Std Dev 0,7013909 Std Dev 0,3164709 
et 4, 
TwrClassRng<127 ||TwrClassRng>=127 TwrClassRng<111 TwrClassRng>=111 
Count 22 || Count 27 Count 24 LogWorth Difference || Count 44 
Mean 18,198485 || Mean 19,201235 Mean 19,291667 18,686089 0,62017 |} Mean 19,746341 
Std Dev 0,6756911 || Std Dev 0,2657808 Std Dev 0,3294337 Std Dev 0,1400155 
PredClassRng<145 || PredClassRng>=145 
Count 7 {| Count 17 
Mean 18,852381 || Mean 19,472549 
Std Dev 0,1214986 || Std Dev 0,1768691 
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G. 


Regression Tree for the Mean 



















































































































































































Number of Smugglers 
Classified 
RSquare N Number of 
Splits 
0,767 128 7 
All Rows 
Count 128 LogWorth Difference 
Mean 1,9338542 21,531724 0,10408 
Std Dev 0,0814334 
TwrClassRng<113 TwrClassRng>=113 
Count 51 LogWorth Difference Count 77 LogWorth Difference 
Mean 1,8712418  12,004726 0,12407 Mean 1,9753247 4,2944776 0,04096 
Std Dev 0,0886991 Std Dev 0,0391438 
PredClassRng<120 PredClass Rng>=120 PredClassRng<90 PredClassRng>=90 
Count 24 LogWorth Difference Count 27 LogWorth Difference | | Count 20 LogWorth Difference || Count 57 
Mean 1,8055556 0,8533217 0,06667 Mean 1,9296296 3,3386588 0,06824 Mean 1,945  1,0583524 0,07778 || Mean 1,9859649 
Std Dev 0,0733377 Std Dev 0,0533761 Std Dev 0,0614398 Std Dev 0,0188429 
TwrClassRng<80 TwrClassRng>=80 PredClassRng<153 |/ PredClassRng>=153 TwrClassRng<139 || TwrClassRng>=139 
Count 12 || Count 12 LogWorth Difference Count 10 || Count 17 Count 5 || Count 15 
Mean 1,7722222 || Mean 1,8388889 1,2311659 0,07048 Mean 1,8866667 || Mean 1,954902 Mean 1,8866667 || Mean 1,9644444 
Std Dev 0,0789237 |) Std Dev 0,0509175 Std Dev 0,0421637 || Std Dev 0,0423956 Std Dev 0,0802773 |] Std Dev 0,04076 
PredClassRng<93_ |/PredClassRng>=93 
Count 7 || Count 5 
Mean 1,8095238 || Mean 1,88 
Std Dev 0,0370899 || Std Dev 0,0380058 
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H. 











Regression Tree for the 





























Mean Number of Smugglers 
Classified with MiniUAV 
RSquare N Number of 
Splits 
0,733 129 6 
All Rows 
Count 129 LogWorth Difference 
Mean 1,9341085 13,06718 0,09795 
Std Dev 0,0893131 
TwrClassRng<110 TwrClassRng>=110 
Count 52 LogWorth Difference Count 77 LogWorth Difference 
Mean 1,875641 21,290582 0,18449 Mean 1,9735931 7,8978124 0,04456 
Std Dev 0,1092887 StdDev  0,038744 
PredClassRng<99 |/PredClassRng>=99 TwrClassRng<143 TwrClassRng>=143 
Count 16 || Count 36 Count 28 LogWorth Difference Count 49 LogWorth Difference 
Mean 1,7479167 || Mean 1,9324074 Mean 1,9452381 8,00203 0,06872 Mean 1,9897959 2,2296143 0,02101 
Std Dev 0,0807087 || Std Dev 0,0619709 StdDev 0,0464495 Std Dev 0,0206426 


























PredClassRng<125 
Count 15 
Mean 1,9133333 
Std Dev 0,0373741 








Count 
Mean 
Std Dev 





PredClassRng>=125 





PredClassRng<96 

















13 Count 13 |} Count 
1,9820513 Mean 1,974359 |} Mean 
0,0220075 Std Dev 0,0308936 


PredClassRng>=96 


1,9953704 
Std Dev 0,0116912 












36 LogWorth Difference 


2,6631592 0,01785 
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TwrHt<27 
Count 
Mean 

Std Dev 





TwrHt>=27 








5 || Count 31 
1,98 |] Mean 1,9978495 
0,0182574 || Std Dev 0,0083244 
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